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INTRODUCTION 


The  present  report  is  the  third  and  last  of  a  series  of 
reports  on  the  theoretical  description  of  polar  motion  and 
nutation.  The  two  previous  reports,  (Moritz,  1980),  henceforth 
referred  to  as  TNP  I,  and  (Moritz,  1981),  referred  to  as  TNP  II, 
will  serve  as  a  basis. 

In  TNP  I,  we  have  considered  the  earth  as  a  rigid  body, 
an  elastic  solid,  and  as  a  "Poincar§  model"  consisting  of  a 
rigid  mantle  and  a  fluid  homogeneous  core.  The  report  TNP  II 
then  treated  the  earth  as  composed  of  an  elastic  mantle  and 
fluid  core,  the  so-called  Molodensky  model. 

The  present  report  considers  in  detail  the  simplest  theory 
of  the  Molodensky  model  due  to  Sasao,  Okubo  and  Saito  (1980). 

An  application  of  Poincare's  equations  of  motion  on  a  Lie 
group,  following  (Moritz,  1982  a), will  provide  a  unified  de- 
ri vati on . 

Finally,  on  the  basis  of  this  model  and  on  the  eigenvalue 
theory  described  in  TNP  I,  secs.  10  and  11,  we  shall  present 
expressions  for  lunisolar  effects  on  polar  motion  and  nutation 
for  the  various  axes:  the  rotation  axis,  the  figure  axis,  the 
"mean  Tisserand  figure  axis",  the  angular  momentum  axis,  and 
the  axis  correspondi ng  to  the  Celestial  Ephemeris  Pole  as  adopt¬ 
ed  by  the  IAU  in  1979. 


THE  SOLUTION  OF  SASAO,  OKUBO  AND  SAITO  (SOS) 


1 .  The  Four  Basic  SOS  Equations 


As  we  already  briefly  mentioned  in  TNP  II1,  pp. 131-133, 
Sasao  et  al.  (1980)  gave  a  particularly  simple  and  elegant 
formulation  of  Molodensky's  problem  in  terms  of  four  complex 
equations,  generalizing  Poincar&'s  (1910)  equations  for  a 
rigid  mantle  and  a  liquid  core. 

The  first  two  equations  may  be  written: 

AO  -  i  (C  -  A)s?  u  +  Ac($  +  inv)  +  Q(£+iiJc)  =  L  ,  (1-1) 

A  ii  +  A  v  +  iC  Qv  +  ni  =0.  (1-2) 

c  c  c  c  x  1 

Here,  A,  A,  C  are  the  (average)  principal  moments  of  inertia 
of  the  earth,  supposed  rotationally  symmetric,  and 


1)  This  notation  is  explained  in  the  introduction:  TNP  i  denotes  (Moritz, 
1980  ),  and  TNP  II  denotes  (Moritz,  1 98 1 ) . 
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Ac,  Ac,  Cc  are  the  corresponding  principal  moments  of  iner¬ 
tia  for  the  core;  fl  denotes  the  (average)  angular  velocity 

of  the  earth's  rotation;  and  the  usual  mathematical  symbols 
2 

are  used:  i  =  -1  and  the  dot  denoting  differentiation  with 
respect  to  time.  The  variables  entering  in  these  equations 
have  the  following  meaning: 


u 

c 


wi  + 


'13 


1u2  , 

+  iC23 


v 


L  -  h  +  ilj 


(1-3) 

(1-4) 


They  are  complex  combinations  of  components  of  the  vectors 


to  =  (uij  ,0)2,013)  ,  £  ~  (xj.X2.X3)  .  (1-5) 

where  to  denotes  the  vector  of  angular  velocity  of  the  rota¬ 
tion  of  the  earth  with  respect  to  inertial  space,  x.  the  vec¬ 
tor  of  angular  velocity  of  the  rotation  of  the  core  with  respect 
to  the  earth's  mantle,  and  l_  denotes  the  lunisolar  torque, 
due  to  the  attraction  of  sun  and  moon.  The  coordinate  system 
(x1  =  x,  x2  =  y,  x^  *  z)  used  is  mantle-fixed  in  the  sense 
that  the  mantle  is  at  rest  (on  the  average)  in  this  system; 
more  precisely,  it  is  a  Tisserand  frame  for  the  mantle  (cf. 

Munk  and  Macdonald,  1960,  p.10;  TNP  1 1,  pp  .  140- 143 ) .  The  x^  axis 
is  directed  to  the  (average)  North  Pole.  Finally,  c^  is  the 
x^^  component  of  the  inertia  tensor  of  the  earth,  and  simi¬ 
larly  for  c2j  ;  c^j  and  c^j  are  the  corresponding  quan¬ 
tities  for  the  core. 

If  the  mantle  is  rigid  and  if  the  coordinate  axes  a  re 
principal  axes  of  inertia,  then  the  inertia  tensor  is  diagonal. 


so  that  the  quantities  (1-4)  are  zero.  Then  (1-1)  and 
(1-2)  reduce  to  Poincarfe's  (1910)  well-known  equations  for 
a  rigid  mantle,  as  they  should.  What  is  surprising,  however, 
is  the  discovery  by  Sasao  et  al.  (1980)  that  the  generalization 
to  Molodensky's  problem,  for  an  elastic  mantle,  is  so  simple. 

The  quantities  (1-3)  and  (1-4)  are  related  by  an 
equation  of  the  form 

c  =  Dn(u-w)  +  012v  ,  (1-6) 

cc  *  D12(u-w)  +  D22v  ,  (1-7) 

where  the  coefficients  ,  012  ,  D22  are  real  constants 

which  depend  only  on  the  elastic  properties  of  the  mantle. 

They  are  easily  expressed  in  terms  of  the  Love  number  k  and 
of  the  well-known  functions  y^(r)  and  y ^ ( r )  in  standard 
notation  introduced  by  Alterman  et  al .  (1959),  or  of  the 
equivalent  functions  R(r)  and  P(r)  used  in.  TNP  II,  Part  B. 
Such  expressions  are  represented  by  eqs.  (54)  to  (57)  of 
(Sasao  et  al .  ,  1980);  they  will  be  derived  in  the  next  section. 
Essential  is  the  fact  that  the  matrix 


D  = 


is  symmetric;  this  is  a  consequence  of  a  reciprocity  theorem 
(Sasao  et  al  .  ,  1980,  eq .  (61)  ). 

The  complex  quantity 


is  related  to  the  lunisolar  (tidal)  potential  Ve  ,  which  we 
assume  to  have  the  form 


V  *  k(xz  cos  ot  +  yzsinot) 

C 


(1-10) 


<  being  a  constant  coefficient.  Such  a  form  correspr  r  to 
a  tesseral  tidal  potential,  being  a  spherical  harmon  of 
degree  2  and  order  1  ,  which  is  responsible  for  i  ced 

nutation  and  polar  motion  (TNP  I,  p.30;  TNP  II,  p.23 
rather  to  the  contribution  of  a  certain  frequency  o  ,  the 
tidal  potential,  which  is  the  sum  of  all  such  contributions. 
Then  w1  and  w2  are  defined  by 


W  =  fi-1«COSCTt 
1  -1 

w2  =  ,7  icsinot 


(1-11) 


which  is  the  desired  relation  between  w  and  Ve  .  Let  it 
also  be  mentioned  that  the  angular  momentum  components  are 
given  by 


L1  =  <MJ2sinot 
L2  *  -<MJ2cosat 


d-12) 


in  view  cf  well-known  relations  (Melchior,  1978,  sec.  2.3; 
TNP  I,  p.30);  here  <  is  the  same  as  in  (1-10),  and 


is  the  well-known  zonal  harmonic  coefficient,  a  denoting 


the  earth’s  semimajor  axis  and  M  its  mass;  C  and  A  are 
(mean)  principal  moments  of  inertia  as  before. 

Eqs.  (1-1),  (1-2),  (1-6),  and  (1-7)  constitute  a  set 
of  4  equations  for  the  4  complex  unknowns  u,v,c,cc  ,  equi¬ 
valent  to  8  equations  for  8  real  unknowns;  they  are  iden¬ 
tical  (apart  from  notation)  to  eqs. (37),  (32),  (54),  (56)  of 
(Sasao  et  al .  ,  1980);  they  will  be  called  the  SOS  equations. 
Now  (1-6)  and  (1-7)  are  simply  linear  algebraic  equations 
with  constant  coefficients,  whereas  (1-1)  and  (1-2)  are 
linear  differential  equations  with  constant  coefficients, 
which  can  be  reduced  to  algebraic  linear  equations  by  seeking 
solutions  of  the  form 

u=uQei0t  ,  v  =  v0eiot  (1-13) 

as  usual.  Thus,  a  linear  system  of  4  ordinary  equations  for 
4  complex  unknown  results. 

This  formulation  of  Molodensky's  problem  is  probably 
the  simplest  given  so  far.  It  is  particularly  remarkable  that 
it  is  valid  for  arbitrary  earth  models  with  an  elastic  mantle, 
a  liquid  core,  and  even  an  elastic  inner  core.  The  detailed 
structure  of  mantle  and  inner  core  enters  only  into  (1-6) 
and  (1-7),  since  the  coefficients  D.^  depend  on  the  elas¬ 
ticity  functions  R (r)  and  P(r)  ;  it  does  not  enter  into 
(1-1)  and  (1  -  2). 

Equations  (1-1),  (1-6),  and  (1-7)  are  relatively 
easy  to  get.  Eq .  (1-1)  is  a  consequence  of  Euler's  equation 
as  generalized  to  a  non-rigid  body  by  Liouville;  it  is  equi¬ 
valent  to  eq .  (10-6)  of  TNP  II.  Equations  (1-6)  and  (1-7) 
will  be  derived  in  the  next  section. 


The  derivation  of  (1-2)  is  considerably  more  difficult. 
Sasao  et  al.  (1980)  obtain  it  by  means  of  the  hydrodynamic 
equations,  which  is  complicated  and  far  from  transparent.  Also 
the  similarity  between  (1-1)  and  (1-2)  remains  unexplained. 
This  is  the  more  regrettable  as  it  is  the  great  merit  of  Sasao 
et  al.  to  have  found  eq.  (1-2),  which  is  simple  and  similar 
tc  (1-1);  an  equivalent  equation  by  Molodensky  (1961;  his 
eq.  (39)  )  is  much  more  complicated. 

For  the  simpler  rigid-mantle  model,  Poincarfe  (1910)  has 
given  two  different  methods  for  obtaining  his  equivalent  to 
our  eq.  (1-2).  One  method  uses  the  hydrodynamic  equations, 
corresponding  to  the  approach  by  Sasao  et  al .  (1980).  Much 
more  interesting,  however,  is  Poincarfe's  other  method,  which 
uses  a  variational  principle  and  exploits  symmetries  expressed 
by  group  theory.  This  method  is  not  only  simpler  and  more 
elegant,  but  also  explains  the  similarity  of  (1-1)  and  (1-2) 
both  are  effects  of  rotating  groups,  the  first  describing  the 
rotation  of  the  earth  with  respect  to  inertial  space,  and 
the  second  expresses  the  rotation  of  the  core  with  respect  to 
the  mantle. 

A  similar  group-theoretic  derivation  of  the  SOS  equa¬ 
tions  will  be  given  in  Part  B  of  the  present  report. 


2 .  Derivation  of  Two  SOS  Equations 


Equati on  ( 1  -  6 ) .  This  equation  describes  the  change  of 
the  inertia  tensor  (elements  c^and  c23)  by  the  deformation 
of  the  earth,  due  to  centrifugal  force,  expressed  by  the 
rotation  vector  u  ,  and  to  lunisolar  tidal  force,  expressed 
by  w  according  to  (1-10)  and  (1-11).  It  is  essentially 
the  same  as  eq.  (63)  of  (Molodensky,  1961)  and  constitutes  a 
generalization  of  well-known  expressions  for  the  rotational 
deformation  of  an  elastic  earth;  cf.  eq.  (14)  of  (Jeffreys, 
1970,  sec.  7.04)  or  TNP  I,  eq .  (4-1). 

We  derive  (1-6)  here  from  eq.(10-32)  of  TNP  II: 


l-i  2 

C13=  -  ^  G  k(<  -  fi  e )C0Sot 
c23=  -  G_1k(<  -  fi2e) sincrt 


(2-1) 


which  can  be  combined  into  one  complex  equation  using  (1-4): 

c  =  -  1  G_1k( <  -fl2e)e iat  .  (2-2) 

Here  G  denotes  the  Newtonian  gravitational  constant  and 
k  the  potential  Love  number  as  usual;  also  the  symbol  a 
for  the  (average)  rotational  velocity  of  the  earth  is  standard 
We  put 


ro  ^  =  Q  to  j  » 


=  ft  ^<*>2 


(2-3) 


so  that  m.  and  in, 


are  the  x  and  x 


componen  ts 


of  the  unit  vector  of  the  instantaneous  rotation  axis  (TNP  I 
sec. 3;  TNP  II,  sec. 10),  and 

m1+im2  =  eeiat  (2  -  4) 

by  eq.  (1-28)  of  TNP  II.  With  (1-3)  this  gives 

u  =  J2ceifft  ,  (2-5) 

and  from  (1  -  11)  we  get 

w  =  a~'<eiat  .  (2-6) 

Thus  (2-2)  may  be  written 

c  =  -  I  G-1k  n(w-u)  .  (2-7) 

3 

In  this  equation  we  have  taken  the  radius  of  the  mean  terres 
trial  sphere,  denoted  by  a  ,  to  be  our  unit  of  length.  If 
a  ^  1  ,  then  we  must  replace  (2  -  7)  by 

c  =  -  ^  G”  1  ka5n( w-u )  ,  (2-8) 

3 

in  order  to  get  the  dimensions  correct.  For  w  *  0  this  re¬ 
duces  to  eq.  (4-1)  of  TNP  I,  as  it  should. 

Following  (Sasao  et  al.,  1980),  we  decompose  the  Love 
number  k  as  foil ows : 


9 


k 


v 

w-u 


(2  -  9) 


with  constant  kQ  and  k^  .  Note  that,  although  u,  v,  w 
are  complex  numbers,  the  quotient  v/(w-u)  is  real  since  the 
common  complex  factor  elot  cancels.  The  constant  k 
represents  the  Love  number  in  the  absence  of  core  rotation 
(v  *  0)  ,  and  k]  expresses  the  effect  of  core  rotation. 

The  decomposition  (2-9)  will  be  justified  later  in  this  sec- 
ti  on . 

The  substitution  of  (2-9)  into  (2-8)  finally  yields 


c  =  Dn(u-w)  +  D12v 

with 

D  =  i  G  " 1  a  5  >1  k  , 

11  3  o’ 

0,2  *  *  is-'a5i)k, 


(2-10) 


(2  -  11) 


This  completes  the  derivation  of  (1-6). 

Equati on  ( 1  -  7 ) .  This  equation  is  the  equivalent  for 
the  core  of  (1-6)  or  (2-10).  Its  derivation,  however,  is 
considerably  more  laborious  and  may  be  skipped  by  the  reader 
who  is  not  interested  in  this  detail. 


Denote  by  C  the  undeformed  core 
o 


FIGURE  2.1.  The  core  before  (CQ)  and  after  defor¬ 
mation  ( c1 ) 

and  by  C1  the  core  after  deformation  (Fig. 2.1  )  .  For  the 

present  purpose  it  is  sufficient  to  consider  CQ  a  sphere 
of  radius  b  ;  this  is  the  usual  spherical  approximation. 

The  xz  product  of  inertia  is  given  by  the  well-known 
formula  (cf.  Heiskanen  and  Moritz,  1967,  p.62): 


The  integral  is  extended  over  the  deformed  kernel  C 

Before  deformation,  this  integral  is  zero  for  reasons 
of  symmetry  (for  the  sphere  as  well  as  for  an  ellipsoid  of 
revolution  if  the  axis  is  the  axis  of  symmetry): 

jpxzdv  =  0  ,  (2-13) 

Co 

where  p  is  the  density  before  deformation,  which  need  not 
be  constant  but  can  be  a  function  of  the  radius  vector  r 
in  the  spherical  approximation: 


p  =  p  (  r ) 


(2  -  14) 


dv  denotes  the  volume  element. 
Since  by  definition 


13 


-Dc 

xz 


(2  -  15) 


( TNP  I,  p . 9 ;  TNP  II,  p.113),  (2-12)  becomes 


-c 


c 

1  3 


(p+Pj)  xzdv 


(2  -  16) 


p1  denoting  the  change  of  density  due  to  deformation.  This 
integral  may  be  split  up  by 


* 


' 

f  f 

= 

+ 

1 

. 

J 

- 

J 

VCo 


as  follows: 


i  3 


pxzdv  + 


I 


pjXZdv  + 


pxzdv  + 


p.xzdv 


n 


c,-c 

1  o 


C  -C 
1  o 


(2-17) 


(2-  18) 


Here  C1-Cq  denotes  the  layer  of  thickness  u  (radial 
component  of  displacement)  by  which  C1  differs  from  CQ 
By  Fig. 2.1  we  have 


C 


1 


u  da 
r 


(2  -  19) 


da  denoting  the  element  of  solid  angle,  or  the  surface 
element  of  the  unit  sphere  a 

The  first  integral  in  (2-18)  is  zero  by  (2-13)  and 
the  last  integral  can  be  neglected  because  it  is  a  second- 
order  quantity,  u  and  p  being  small  of  first  order. 


Thus  there  remains 


c 

C  13s 


p1  xzdv 


b2jjp  xzurda  ; 


(2  -  20) 


here  we  have  used  (2-19). 

By  the  well-known  Poisson  equation  we  express  the  ano 
malous  density  p^  in  terms  of  the  corresponding  potential 
disturbance  V1  : 


-4irGp  1  ; 


(2  -  21) 


cf.  TNP  II,  eq.  (6-30)  .  Thus  the  first  integral  of  (2-20) 
may  be  written 


p.xzdv  » 


i_(f[ 

[7TG  Jjj 


aV1  xzdv 


(2  -  22) 


This  integral  can  be  transformed  by  means  of  Green’s 
second  identity: 


UAVdv 


[  [ V a U d v  +  f  j  'u  ^  -  V  IdS  ;(2  -  23) 
> '  JJlan  an  i 


cf.  (Heiskanen  and  Moritz,  1967,  p.ll).  Here  v  denotes 
a  volume  bounded  by  a  closed  surface  S  ,  and  3/an  is  a 
symbol  for  the  derivative  along  the  surface  normal.  In  the 
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present  case,  v  is  the  spherical  core  Cc  ,  the  boundary 
S  is  the  sphere  r  *  b  ,  the  normal  derivative  is  taken 
along  the  radius  vector  r  : 


I _  ,  3_ 

3  n  3  r 

and  the  surface  element  is 
dS  *  b2d* 


(2-24) 


(2  -  25) 


by  Fig.  2.1;  furthermore  we  take  in  (2-23) 


V  *  Vt 


(2-26) 


U  *  xz  =  rzsinecoseco$> 


(2-27) 


in  spherical  coordinates  r,  e,  x 
Thus  (2  -  23)  gives 

III  ■  b1j(*z  TF1  ■  v, 

C  a 

o 

the  first  integral  on  the  right-hand  side  is  zero  since 
U  a  xz  is  a  harmonic  function  (AU*0) 

The  differentiation  of  (2-27)  shows  that 


(2  -  28) 


=  2rsi  n  0CO-9  ;:os  *  *  ^r- 


Hence  by  (2-22)  and  ( 2 - 2C )  we  may  transform  (2-2C)  into 


b2  (i  i  3Vi  2V,  ) 

13  3  j  J  xz  !  lr~  "  ~r~  “  ^Gpu?j  ^  * 


(2  -  29) 


For  the  final  transformation  we  must  now  use  the  for¬ 
mulas  for  spherical  elasticity  derived  in  TNP  II,  Part  B; 
they  may  also  be  applied  to  a  liquid  core  (homogeneous  or 
with  a  density  law  (2 -14))  since  the  equations  of  hydro¬ 
dynamics  may  formally  be  considered  a  special  case  of  the 
equations  of  elasticity  (TNP  II,  p.63). 

By  TNP  II,  pp.82,  85,  and  94  we  put 

ur  =  H  (  r )  S(  e  ,  A  )  ,  (2-30) 

Ve  +  v,  R(  r)S(  8,  A)  ,  (2-31) 

where  S(e,A)  is  defined  as  the  value  of  the  lunisolar 
potential  Vg  at  the  earth's  surface  (r  =  a)  : 

S  ( 6 ,  a  )  =  Ve(a,9,A)  •  (2-32) 

We  wri te  ( 1  -  10 )  as 

=  <(xzcosat  +  yzsinat)  = 

(2-33) 

2 

=  <r  sin9cos8(cosAcosot  +  sinAsinot)  , 
o  r 

1  2  “  "> 

V0  =  -  <r  R2 1 ( e ,A )cosat  +  S21 (9  ,  A  )  s i no t  j  ,  (2-34) 

using  Legendre  surface  harmonics  in  the  notation  of 
(Heiskanen  and  Moritz,  1967,  p.29). 


From  (2  -  34)  we  get  immediately 


3 Ve  _  2 Ve 
3r  r 


=  0 


so  that  by  ( 2  -  31) 


3  V  i  .  2V_ 
3  r  r 


[R,(r)-2r’1R(r)]  S(e.A)  ,  (2  -  35) 


putting  R'(r)  =  dR/dr  .  Thus  the  expression  between  pa¬ 

rentheses  in  (2-29)  becomes 


3V,  2V_. 

3r  r 


[r'2P(r)' 


-  4irGpuf  =  p  R,(r)-2r”1R(r)-4irGpH(r)J  S  (  0  ,  X  ) 

(2  -  36) 


2r"1R(r)]  S(0  ,X)  , 


introducing  the  function  P ( r )  defined  in  TNP  II,  p.88, 
eq.  (7-40). 

We  now  put 


xz  =  -  r2  R  (0,x) 

3  1 1 


(2  -  37) 


(cf.(2-33)  and  (2-34)  )  and  note  that  the  integral  (2-29) 
is  extended  over  the  sphere  r  *  b  .  Substituting  (2-36) 
and  (2-37),  with  r  =  b  ,  into  (2-29)  we  have 


(2  -  38) 


P(b)-2bR(b)  I 
1 2  tr  G  L  J  >> 

a 


R21 (0 ,A)S(e  ,X)do  . 


By  (2-32)  we  get  $(e,A)  by  putting  r  =  a  in  (2-34). 
We  have 


(Heiskanen  and  Moritz,  1967,  p . 2  9 ) .  Thus 


R21(e,x)S(e,x)do  * 


4  it  .2  4. 

—  a  < t o So t 


and  (2-38)  takes  the  final  form 


■c  .  aV 
'13  15G 


|p(b)-2bR(b)J  < cosa t 


(2  -39  ) 


The  quantity  is  obviously  given  by  the  same  express¬ 

ion,  with  cosat  replaced  by  sinot 


In  the  expression  (2-33), 


denotes  the  potential 


of  al 1  perturbing  forces  acting  on  the  body.  This  includes 
the  lunisolar  tidal  potential  and  the  disturbing  centrifugal 
potential.  If  we  want  to  restrict  k  to  the  effect  of  the 
lunisolar  potential  only,  we  roust  add  the  effect  of  the 
centrifugal  perturbation  (TNP  II,  p.115),  replacing  < 


by  K-n  e 


Thus  from  (2-39)  we  have 


+  c. 


aV 


jp(b)-2bR(b)J  (ic-fl2e  )e  ' ' 


(2  -  40) 


or,  by  (2-5)  and  (2-6),  finally 


aV 


jp(b)-2bR(bj]  ft(w-u)  . 


(2-41) 


This  formula  expresses  the  elements  and 

of  the  inertia  tensors  of  the  core  deformed  by  the  effect  of 
lunisolar  (w)  and  centrifugal  (u)  perturbati on . 

As  a  check  we  replace  b  by  a  ;  then  the  core  is 
replaced  by- the  whole  earth,  and  (2-41)  should  reduce  to 
(2  -  7).  In  fact,  putting  a  =  b  =  1  in  (2  -  41)  we  have 


P ( a ) -2R( a )  a(w-u) 


(2  -  42) 
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By  TNP  II,  eqs.  (8-20)  and  (8-24),  we  have 


P ( a )  =  2  -  3  k  ,  R( a ) 


1  +  k  , 


(2  -  43) 


which  on  substitution  into  (2-42)  gives  (2-7)  as  it  should 
be . 

As  a  final  step  we  perform  a  decomposition  analogous 
to  (2-9)  : 


P(b)  =  PQ(b)  +  -i-  P,  (b)  , 


R(b)  =  R0(b)  +  R -j  ( b ) 

w-u 


(2  -  44) 


The  possibility  of  this  decomposition  is  a  consequence  of 
the  linear  character  of  the  partial  differential  equations 
of  elasticity  and  of  the  boundary  conditions.  In  fact,  by 
eqs.  (8-25)  of  TNP  II,  the  elastic  functions  P(r)  and 
R(r)  depend  linearly  on  the  Love  numbers  h,k,l  , 
which  by  (10-44)  of  TNP  II  are  the  solution  of  three  linear 
equations  of  the  form 


a.h  +  a.k  +  a,l  =  a^  , 
1  Z  3  o 


bh  +  bk  +  bl=b  +  b,- 

1  2  3  o  4, 


c.h  +  c.k  +  c . 1  =  c 
1  2  3  o 


(2  -  45) 
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In  (10-44)  of  TNP  II  we  have  written  the  second  equation 
in  a  slightly  different  form,  making  them  dependent  on  the 
Molodensky  parameter  e  .  However,  as  eqs.(lQ-43),  ibid.  , 
show,  the  dependence  is  rather  on  .  If  we  have  both 

tidal  and  centrifugal  perturbati on ,  we  must  replace  this 
factor  by 

— 1 T—  ■  (2-«) 

K  “Q  £. 

By  eq.(4-60)  TNP  II,  8  is  proportional  to  our  v  , 
(this  also  holds  for  an  elastic  mantle),  so  that  the  quotient 
(2-46)  is  proportional  to 

(2-47) 

w-u 

(a  real  quantity,  since  the  complex  factor  eICTt  cancels!), 
which  proves  (2-45).  Of  course,  the  coefficient  b^  is 
now  different,  but  all  other  coefficients  are  equal  in 
(2  -  45)  and  eq. (10  -  44)  of  TNP  II  . 

The  solution  of  (2-45)  now  does  give  h  ,  k  ,  1  as  linear 
functions  of  the  ratio  (2-47),  of  the  form  (2-9),  where¬ 
upon  eqs .  (8-25)  of  TNP  II  give  P(b)  and  R(b)  in  the 
form  (  2  -  44 ) . 

The  substitution  of  (2-44)  into  (2-41)  thus  gives 

cc  =  D2 , ( u-w)  +  D22v  (2-43) 


where 


°21  =  -  ^  a2b2G-1«[po(b)-2bRo(b)] 

d22  -  —  a2b26“1fl[pi(b)-2bR1(b)] 


(2  -  49) 


Equations  (2  -  48)  and  (2  -  49)  are  equivalent  to  equations  (56) 
and  (57)  of  (Sasao  et  al.,  1980). 

We  finally  mention  the  symmetry  property 


(2  -  50) 


Comparing  the  corresponding  expressions  given  by  (2-11) 
and  (2  -  49),  the  equality  (2  -  50)  is  far  from  evident.  In 
fact,  this  equality  is  a  consequence  of  a  deep  theorem  of 
elasticity  (Betti's  reciprocity  theorem).  We  shall  not  prove 
it  here,  referring  the  reader  to  (Sasao  et  al  .  ,  1980,  sec. 5). 

Eq.  (2-48),  with  the  symmetry  (2-  50).  is  indeed  equi¬ 
valent  to  (1-7). 
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PART  B 

APPLICATION  OF  POINCARE'S  VARIATIONAL  PRINCIPLE 


3 .  The  Rotation  Group 

The  elementary  properties  of  the  rotation  group  will 
play  a  basic  role  in  Poincar&'s  variational  principle,  as 
already  indicated  at  the  end  of  sec.l.  They  will  therefore 
be  dicussed  in  this  section,  following  (Moritz,  1982b). 

Let  a  rotation  in  R^  (threedimensional  Euclidian 
space)  be  represented  by 

x'  =  Ax  ,  (3-1) 

x  and  x'  being  vectors  and  A  denoting  a  3x3  matrix. 
Rotation  matrices  have  the  properties 

1.  The  product  of  two  rotation  matrices  A  and  13  is  again 
a  rotation  matrix  C  =  AB 

2.  For  the  unit  matrix  _I_  we  have 

A  J.  *  I A  »  A  .  (3-2) 


3.  Every  rotation  matrix  A  has  an  inverse  A”  which 
again  is  a  rotation  matrix. 

This  properties  characterize  the  mathematical  structure  of  a 
group ;  we  therefore  speak  of  the  rotation  group. 

Another  important  property  of  rotation  matrices  A 
is  that  the  inverse  is  simply  the  transposed  matrix: 

A"1  =  AT  ,  (3-3) 

so  that  we  have 

A  A  T  =  A  T  A  =  l  .  (3-4) 

Denote  the  elements  of  A  by  a..  : 

i  j 

A  »  [  a  ]  .  (3-5) 

Any  rotation  can  be  fully  described  by  means  of  the  three 

variables  q  ,  for  which  we  may  take  the  Euler  angles  * , 

r 

e,  ip  (TNP  I,  p.51).  Thus  the  a  will  be  functions  of 

'  j 

qr  : 

a.  .  »  a.  .  (q  ,q  ,q  )  =  a.  .  (q  )  .  (3-6) 

i  j  i  J  1  2  3  i  j  r 

These  functions  are  easily  seen  to  be  continuous  and 
differentiable.  Thus  the  rotation  group  is  a  continuous  group, 
or  Lie  group. 

An  infinitesimal  change  of  q  changes  the  matrix  A 


r 


*4 


3  A 

dA=  — dq  (3-7) 

3<r 

(summation  convention!). 

The  matrix 

dn  =  A"’  d*  "  ATdA  (3-8) 

will  be  skew-symmetric,  which  is  an  immediate  consequence  of 
differentiating  (3-4)  : 

ATdA  +  dATA  *  ATdA  +  (ATdA)  T  *  0 

Thus  it  has  the  form 


1 

0 

-d*3 

d*2 

* 

1 

i 

d*3 

0 

_dir1 

(3-9) 

1 

1 

1 _ 

-d,2 

dwl 

0 

_ 

j  On  introducing  the  matrices 

(3-10) 


I  0 

1 

0 

o' 

r 0 

0 

1 

ro  -1 

1 

o'! 

1 

4 

£i  - !  0 

0 

-1 

,  E2  =  0 

L  \ 

0 

0  1  *  £3 

1 

»  1  1  0 

1 

0  !  , 

i 

Lc 

1 

0. 

L-i 

0 

OJ 

0  0 

1 

this  may  be  written 


dn  =  d*  +  £2  d*2  +  £3  dlT3  *  £i  dlT  j  •  ( 3  -  11 ) 

The  geometrical  interpretation  of  the  term  Ejdir,  is  clear 
it  represents  a  rotation  by  the  infinitesimal  angle  d^ 
around  the  x1  axis,  and  similar  for  the  other  terms. 

As  (3-  11)  shows,  infinitesimal  rotations  are  commu¬ 
tative  (finite  rotations  are  not:  there  is  AJ3/*  BA  in 
general  ) . 

The  angular  velocity  component  w1  may  be  considered 
a  change  of  d* 1  with  respect  to  time: 


01 


1 


(3-12) 


and  similarly  for  u2  and  u»j  .  Denoting  dn/dt  by 
n  ,  we  have  from  ( 3  -  11 ) : 

n  *  Eju,  .  (3-13) 

The  matrices  £.  satisfy  the  basic  commutation 
relations 


t£, .  £2^  =  £3 

t£2*  h]  -  £1 
£3  *  £p  -  £2 


(3  -  14) 


Here  the  commutation  symbol  [}  stands  for 
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Hi  -  £j] 


E.E.  -  E.E. 
-■-J  -J-i 


(3-  15) 


being  the  usual  matrix  product  of  E.  and  £. 

The  matrix  dii  is  skew-symmetric  and  directly  express¬ 
ed  in  terms  of  three  infinitesimal  parameters  dir.  .  It  is 
thus  simpler  than  dA  .  The  latter  can  be  obtained  by  (3-8), 


dA  »  Adn 


(3-16) 


which  expresses  dA  also  in  terms  of  the  basic  "infini¬ 
tesimal  group  variables"  ditj  ,  d-n 2  ,  du^ 

It  is  clear  that  the  matrix 


l  +  dn 

represents  a  small  rotation  which  is  close  to  the  unit  matrix. 
Thus  we  may  interpret  (3  -  S)  by  saying  that  the  multiplication 
by  A-1  =  AT  transforms  an  arbitrary  small  rotation 

dA  into  a  matrix  dn  "in  the  neighbourhood  of  the  unit 

matrix". 

For  a  general  Lie  group,  the  commutation  relations 
(3-  14)  are  replaced  by 


f£« 


E  . 
-J 


=  C.  ..  E, 
i  j  k— k 


(3  -  17) 


where  the  c,-jk  are  constant5»  called  the  structure  constants 
of  the  group.  For  a  general  group  the  indices  i,  j,  k  run 
from  1  to  n  ,  n  being  again  the  number  of  degrees 
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of  freedom.  By  (3-  15),  the  interchange  of  i  and  j  means 
a  change  of  sign  .  Hence, 


cjik  3 

“C  i  jk  ' 

(3  - 

18) 

the  rotation  group  we  have  in  particular 

C  1  2  3  = 

C  2  3  1  =  C  3  1  2  =  1 

C  2 1  3 

C 32 1  *  C 1 32  =  _1 

(3- 

19) 

al  1  other  c .  =0  . 

i  jk 


Invariant  differential  forms. By  substituting  (3-7)  into 
(3-8)  we  find 

3  A 

dn  =  AT — -  dq  .  (3-20) 

3qr 

In  view  of  (3-9)  this  has  the  form 

dtr.  =  a  r .  dq  r  ,  (3  -  21) 

expressing  dir.  as  a  linear  combination  of  dq  ;  the 
summation  convention  holds  as  usual. 

If  we  replace  A  by  E3A  ,  with  an  arbitrary  con¬ 
stant  rotation  matrix  B  ,  the  right-hand  side  of  (3-20) 
becomes 


< 
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T  3A  _  T  3A 

(B  A  )  B  —  dq  «  a'b'b  —  . 

3qr  aqr 

identical  to  (3-20)  since  BTB  =  _I_  for  a  rotation  £  . 

Thus  a  multiplication  of  all  matrices  A  by  the  same  con¬ 
stant  rotation  matrix  B  from  the  left  ("left  translation 
of  A  by  B  ")  leaves  the  form  (3-20)  or  (3-21)  invariant. 
Hence  we  speak  of  an  invariant  (more  precisely,  left-invariant) 
differential  form  for  the  group  under  consideration. 

The  coefficients  ari  are  functions  of  qr  .  Note 
particularly  that,  in  general,  the  differentials  dir.  cannot 
be  integrated  to  give  new  coordinates  * .  .  Mathematically 

speaking,  the  d  tt  .  are  not,  in  general,  "perfect  differentials". 
One  also  speaks  of  nonholonomic  coordinates,  which  make  sense 
only  in  the  infinitesimal  domain;  cf . ( Graf arend ,  1975). 

Briefly,  the  qr  are  holonomic  but  not  group  inva¬ 
riant,  whereas  the  diT.  are  group-invariant  but  not  holonomic. 
The  property  of  being  group-invariant  is  so  important,  however, 
that  the  du.  are  basic  in  the  theory  of  continuous  groups. 

Dividing  by  dt  and  noting  (3-12)  we  get  from  (3-21) 


The  inverse  relation  may  be  found  by  solving  (3-22)  for  qr  . 
It  has  the  form 


t-v; 

>  -j 
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■ 

^  r  3  Sri“i 

(3  -  23) 

or,  corresponding  to  (3-21), 

m 

dq  ~  $  .dir. 

^  r  wri  i 

(3-24) 

i 

i 

i  / 

It  is  clear  that  the  matrix 

rs  .] 

r  i  J 

is  the  (transposed)  in- 

verse  of  the  matrix  [a  .] 

r  i 

,  so  that  the  relations 

s 

a  .  8  .  *  6  .  . 

n  rj  ij 

a  .  8  . 
MSI 

=  5rs  (3  -  25) 

V 

hold;  the  Kronecker  deltas 

of  the  unit  matrix. 

6 .  .  and 
'  J 

6rs  denote  the  elements 

To  fill  these  abstract  formulas  with  a  concrete  meaning. 

1 

let  us  note  that  (3-22)  is 

simply  a  general  form  of  Euler’s 

■MB 

kinematical  equations  (TNP 

I,  P.51) 

so  that  for  the  rotation 

group  the  matrix  [ cx  r ;  ] 

has  the  form 

»_ 

]-  -  s  i  n  e  s  i  n*  - 

cos*  0 

• 

.  ,  i 

[o  r . j  =  -s i necos  * 

sin*  0 

>  (3-26) 

|_  cose 

0  1  _j 

taking  q  ?  *  -j,  ,  q 2  *  9  , 

q,  =  * 

.  The  explicit  computa- 

• 

tiort  of  the  matri  x  r  8  .  , 

r  i 

by  inverting  (3-26)  is  left  as 

h 

an  exercise  to  the  reader. 

- 

The  generators  of  a  Lie  group 

Let  us  form  the  differ- 

■;  - 

ential  of  a  function 

• 

► 

> 

■ 

• 

*.  -  -  -  - 

-  - * 

(3*27) 


(3  -  28) 


3f 


(3-29) 


Note,  however,  that  the  notation  df/jiTj  has  only  a  formal 
character,  it  is  not  a  true  partial  derivative  of  f  with 
respect  to  n .  since  for  nonholonomic  coordinates  only  dn . 
but  not  tt  j  exists. 


In  Lie  group  theory  it  is  customary  to  write 


instead  of  af/3ir 


so  that  symbolically. 


.  ®  r  i  „  _ 

3-n- .  aq  r 


(3  -  30) 


by  (3-29).  The  differential  operators  X.(i  =  1,  2,  3)  are 
called  the  generators  of  the  group. 

The  torque.  If  we  consider  the  rotation  group,  and  if 
the  function  f  is  the  potential  energy  U  ,  then  the 

generators  X.  have  an  important  physical  meaning.  In  fact. 


L.  =  X.U  = 


(3  -  31) 


is  nothing  else  than  the 
Let  us  consider  the  work 


component  of  the  torque  L 
done  by  a  small  rotation 


about  the 


axis.  We  have 


=  L 1  d  ir  , 


(3-32) 


this  expression  is  the  rotational  analogue  for  the  work  done 
by  a  small  translation  d x 1  along  the  x ]  axis: 


dW  =  K  dx. 

t  ran s  1  1 


K  being  the  force  component  along  the  x^  axis.  The 

change  of  potential  energy  dll  is  equal  to  the  work  dW 
Taking  into  account  also  small  rotations  about  the  other  co- 


ordinates  axes 


,  ( 3  -  32 )  gives 


x„  and  x, 

2  3 

dU  =  L  J  dir  J  +  l^diTj  +  L  3  dir  3  =  L^dir^ 

On  the  other  hand,  (3-28)  and  (3-30)  yield 


dU  =  -iiL  dir,  =  X,  Udrr,  • 

,  k  k  k 

3irk 

The  comparison  of  these  two  expressions  gives  (3-31),  which 
was  to  be  shown. 

Commutation  relations  and  structure  constants  again. 
Consider  the  differential  of  a  rotation  matrix  A  .  By  (3-28), 
(3-30),  (3-16),  and  (3-11)  we  have 

dA  =  X.Adi,  .  =  A  £  .dir.  • 

-  i  —  i  —  i  i 


wnence 


X . A  =  A  E  .  (3-33) 

This  means  that  applying  the  differential  operator  X  to 
a  rotation  matrix  A  is  equivalent  to  multiplying  A 
by  the  matrix  E. 

—  i 

A  s trai ghtforward  consequence  o*  this  fact  (please  veri¬ 
fy!)  is  the  basic  theorem  that  the  generators  X.  satisfy 
the  same  commutation  relations  (3-17)  as  the  matrices  E.  : 


X.  ,  X.  =  C.  ..X,  , 

i  j  i  j  k  k 


(3  -  34) 


where 


[X.  ,  X.l  =  X.X.  -  X.X.  . 
i  J '  '  J  j  ■ 


This  must  be  interpreted  according  to  (3-30) 


[X.  ,  X  ]f  =  X.X.f  -  X.X.f 
'  J  '  J  J  ' 


(3-35) 


=  _L_  fg  ±L  i  -  8  -L-  fa  — 

riaqr  [  s  J  j  rj3qr  f  siaq, 


The  differentiation  gives 


rjsqr  (  s  i  3 q s 


(3-36) 


[X.  ,  X.  ]f  =  j  s  .-5^- 
'  J  !  r  i  3q  r 


36  .  'i  .f 
s  i  j  _3j_ 

rj  r  J  sq 


since  the  second  derivatives  3  </3q  3qs  cancel.  In  view 
of  (3-25)  this  may  be  written 


X. ,  X  . 1 f  = 


.  k  .  r  i  <  q 


i!»i  }o  ii- 

r  j  3 q  r  J“pk3qp 


The  last  two  factors  are  X  f  ,  by  (3  -  30).  Hence  the  compa 


rison  with  (3-34)  gives  the  structure  constants: 


c .  . .  =  a  , 

i  j  k  sk 


36 


11 


1.  r  '  3q  r 


(3  -  37) 


expressed  in  terms  of  the  functions  a  .  and  e  . 

r  i  r  i 

An  alternative  expression  is  obtained  by  differentia¬ 
ting  the  first  equation  of  (3-25)  written  in  the  form 


a  .6  .  =  6  . . 

sk  s  j  j  k 


obtaining 


36 


sk3q 


11 


3a 


+  6 


sk 


sj  3q  r 


=  0 


By  means  of  this  relation,  (3-37)  is  easily  transformed  into 


da 


Cijk  =  Sri6sjl3q 


rk 


3a  . 

W7) 


(3  -  38) 


As  an  exercise,  the  reader  is  invited  to  compute  the 
structure  constants  (3  -  19)  of  the  rotation  group  from  (3  -  38), 
using  (3-26). 

It  is  surprising  that  the  result  of  evaluating  (3-38), 
using  the  functions  i  .  and  3  .  ,  is  a  constant.  In  fact, 

nonhol onomi c  coordinates  may  be  introduced  by  an  equation  of 
form  (3-21),  with  a  .  being  arOitrary  functions  of  q  , 


also  for  general  dynamical  systems  not  related  to  Lie  groups. 
Then,  of  course,  the  cijk  computed  by  (3-38)  wi  1 1  not  in 
general  be  constants. 

The  fact  that  c...  are  constants  if  the  dir.  are 

i  j  k  i 

invariant  forms  for  a  Lie  group  reflects  an  essential  symme¬ 
try  of  the  group:  the  relation  (3-33)  is  the  same  at  every 
point  of  the  group  space. 

Let  it  finally  be  mentioned  that  the  structure  constants 
are  all  zero  if  the  group  is  commutative  (Abelian);  then,  by 
(3-38)  we  have 


(3-39) 


These  are  the  integrabi  1  i  ty  conditions  for  dir.  ,  which 
can  then  be  integrated  to  give  true  (holonomic)  coordinates 

jr  . 

\ 

More  about  Lie  groups  can  be  found,  e.g.,  in  (Smirnow, 
1971)  or  ( Choquet-Bruhat  et  al.,  1977). 


According  to  Hamilton's  principle,  well  known  from  clas- 
sical  mechanics,  a  dynamical  system  moves  from  time  tQ 


to  time 


1 


in  such  way  that 


( T  -  U)dt  =  extremum 


(4-  1) 


T  being  the  kinetic  energy  and  U  the  potential  energy. 
Here  U  is  assumed  to  depend  on  n  parameters  (gener¬ 
alized  coordinates)  q.  (i  =  1,2, . ,n),  whereas  T 

depends  on  q.  and  on  the  generalized  velocities 


q:  = 


dq; 

W 


If  we  introduce  the  lagrangian  function 

£  =  T  -  U  (4-2) 

which  is  a  function  of  q.  and  q.  ,  then  the  variational 

i  i 

principle  (4-1)  leads  to  Lagrange's  equations 


All  this  can  be  found  in  any  text  on  classical  mechanics,  such 
as  (  Goldstein,  1980),  (Lanczos,  1970),  or  (Whittaker,  1961). 

Now  following  Poi nca  r§  (1901),  we  shall  introduce  non- 
holonomic  "velocities"  to.  by  (3-22),  considering  £  a 


function  of 


and 


0)  . 


£  =  £  ( q  r ,  <o  .  ) 


(4-4) 


To  make  a  distinction,  we  shall  denote  £  as  a  function  of 


q  • 


and  g .  by  £ 


£=  £(qr>  to.)  =  £(qr,  i\  r) 


(4-  5) 


Thus  the  Lagrangian  equations  (4-3)  must  now  be  written 


d 

oT 


3  £ 


=  0 


(4  -  6) 


as  we  did  in  the  preceding  section,  we  shall  use  the  subscripts 
i,  j,  k,  ...  for  the  du's  and  the  to '  s  and  r,  s, 
for  the  q's  and  g's  ;  both  sets  of  subscripts 

run  from  1  to  n  ,  n  being  called  the  number  of  degrees 
of  freedom  for  the  motion. 

The  differentiation  of  (4-5),  using  (3-22),  gives 


so  that 


d 

=  arkdT 


+ 


3£_  "“rk 
3uk  dt 


(4-7) 


Now  we  multiply  (4-6)  by  ,  defined  by  (3-23),  and 

substitute  (4-7).  Noting  that 


S.o,  ■  6  .  , 
r  i  rk  i  k 


by  ( 3  -  25 ) ,  we  obtain 


Now 


and 


*  UsJ  + 

at  [aw . J  "r  i 


dark 

~TT 


sq. 


L. 


_  3  £ 


3d 


3Q 


1£_  ^rk  ,  1£_ 

k  -  rl9<r 


6':3w.  “3T 


9a  rk  dqs 


3  £ 


3a 


di¬ 


rk 


3a 


sq, 


=  3 


rk 


sj3q, 


“0)  . 
J 


3  w , 


3  ai,  3d 

k  r 


jot 


s  k 


(4-8) 


(4  -  9) 


(4  -  10) 
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by  (3-22).  The  substitution  of  (4-9)  and  (4-10)  into  (4-8) 
gives 


d  i  3£  ; 

■JT  3  a) .  J 


+  6ri  [  3  q 


r  i  3  q 


(4  -  11) 


Finally  we  substitute  (\ 


from  (3-23) 


^  '  8Sj“j  ' 


obta i ni ng 


d  f d £  x  0  fo0t 
hi  +  6 r i  6 s j  [Jq, 


3Ctsk)  3  £ 

3qr  j  “j3«k 


*rlltr  *  °  -(4-12) 


If  the  anholonomic  velocities  u .  and  the  corresponding 
anholonomic  "coordinates"  dw.  are  group  variables  in  the 
sense  of  the  preceding  section,  then  we  may  introduce  the 
structure  constants  of  the  group  by  (3-38)  and  the  generators 
X.  =  3  /  3  Tf .  by  (3  -  30).  Then  (4  -  12)  reduces  to 


d  |  3  £ 
d  t  3  ui . 


.  -  3  £  3  £  _  n 

+  c  .  . .  oj  .  — -  -  — -  =  U 

i  j  k  j  3uk  3 tt  . 


(4  -  13) 


These  are  Poincare's  (1901)  equations  of  motion  on  a 
Lie  group.  We  shall  also  speak  of  Poincare's  variational  prin¬ 
ciple  although  it  is  a  new  formulation  of  the  classical  prin¬ 
ciple  (4-1)  rather  than  a  new  variational  principle.  It  will 
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be  seen  to  be  fundamental  for  all  questions  of  earth  rotation 
involving  a  liquid  core;  in  fact,  Poincare's  mathematical  in¬ 
vestigations  were  motivated  precisely  by  the  liquid-core  pro¬ 
blem. 

Most  classical  textbooks  on  theoretical  mechanics  do  not 
give  it.  An  exception  is  (Whittaker,  1961,  pp. 42-43),  who  con¬ 
siders  general  anholonomic  coordinates  (not  necessarily  on  a 
group,  so  that  the  c_k  are  not  necessarily  constants) 
and  whose  derivation  we  have  been  following. 

Nowadays,  thi s  topic,  motion  on  a  Lie  group,  is  quite 
fashionable,  mainly  due  to  the  work  of  Arnold  (1978,  Appen¬ 
dix  2);  see  also  (Hermann,  1968,  ch.16  and  33)  and  (Abraham 
and  Marsden,  1978,  ch.4).  All  these  treatments  have  a  consi¬ 
derable  level  of  sophistication;  most  accessible  is  Hermann, 
who  on  pp. 171-172  gives  a  modern  derivation  of  (4-13)  in  terms 
of  external  differential  forms,  and  most  difficult  is  Abraham. 
None  of  these  authors,  not  even  Whittaker,  however,  mentions 
Poi ncarS ! 

We  note  that  if  the  structure  constants  c...  are  zero, 

i  j  k 

then  the  i  ntegrabi  1  i  ty  conditions  (3-3y)  are  satisfied:  we 
nave  true  (hoionomic)  coordinates  * .  with  u.  =  A.  .  Then 
Poincare's  equations  (4-13)  reduce  to  the  Lagrangian  equa¬ 
tions  (4-3),  with  it.  instead  of  q.  ,  as  it  should  be. 

To  return  to  the  general  case,  we  see  that  Poincare's 
equations  (4-13)  differ  from  the  Lagrangian  equations  (4-3) 
only  by  the  second  term  in  (4-13),  which  expresses  the  non- 
holonomity  of  the  group  variables  and  involves  the  structure 
constants  of  the  group. 

Taking  into  account  (4-2)  and  the  fact  that  U  does 
not  depend  on  qr  and  hence  is  independent  of  *i.  ,  we  may 


write  (4-  13)  in  the  form 


C  .  .  .  to 
i  jk 


3  T 

j3o,k 


(4  -  14) 


we  are  from  now  on  using  the  opposite  sign  convention  for  L.  : 

-L.  =  =  X.U  (4  -  15) 

by  (3-31).  This  is  the  form  derived  by  Poi ncar§(  1901 )  and 
used  by  him  (PoincarS,  1910)  in  his  elegant  treatment  of  the 
liquid-core  problem  for  a  rigid  mantle  (his  c;jk  have  diffe¬ 
rent  sign). 

Application  to  Euler's  equations.  Let  us  consider  the 
rotation  of  a  rigid  body.  Then  the  group  under  consider¬ 
ation  is  the  rotation  group  in  ,  whose  structure  constants 

are  given  by  (3-19).  The  kinetic  energy,  for  principal  axes 
of  inertia,  is  well  known  to  be 

T  =  ( Am  i 2  +  Bu)22  +  Cu)32)  ,  (4  -  16) 

as  any  textbook  on  analytical  mechanics  shows.  Then  (4-14) 
gives  immediately 


A  i  ^ 

+ 

(C  - 

c<\ 

3 

CM 

3 

CO 

=  L]  ’ 

8  uj  j 

+ 

(A- 

C )  o>  2  U)  j 

=  L  2  ’ 

C  j . 

+ 

(8  - 

A  )  W  -  ul  - 

=  L7  ’ 
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that  is,  Euler's  equations  for  the  rotation  of  a  rigid  body 
(cf.  TNP  I,  p.10). 

The  following  sections  will  apply  (4-14)  to  more  rea¬ 
listic  earth  models:  the  Poincarb  model  (rigid  mantle  and 
liquid  core),  and  the  Molodensky  model  (elastic  mantle  and 
liquid  core),  leading  us  to  the  theory  of  (Sasao  et  al . , 
1980) . 


5. 


This  model,  the  Poincar§  model,  has  been  treated  in 
TNP  I  ,  secs. 12  and  13.  We  shall  here  derive  the  basic  equa¬ 
tions  (12-8)  and  (12-9)  of  TNP  I,  which  were  given  there 
without  proof. 

We  recall  some  basic  mathematical  features  of  the 
Poincar&  model,  following  sec. 12  of  TNP  I 

Let  us  refer  the  ellipsoidal  shell  to  principal  axes 
xyz  ;  then  the  inner  ellipsoidal  surface,  which  encloses  a 
liquid-filled  cavity,  has  the  equation 


By  the  change  of  variables 


x  1 


(5  -  2) 


this  surface  is  transformed  into  the  unit  sphere 

x|2  +  y'2  +  z'2  =  i  .  (5-3) 

Poincare  considers  a  motion  of  the  liquid  such  that,  by  the 
transformati on  (5-2),  it  is  transformed  into  a  rotation  of 
the  sphere  (5-3).  Thus  the  velocity  in  the  auxiliary  x'y’z' 
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system  is 


X2Z  ' 

-  x3y‘~ 

^  X  X  = 

V  1 

x3X 

-  xiz‘ 

xly' 

X 

Csl 

X 

f 

(5-4) 


if  the  corresponding  rotation  vector  is  ^  .  Going  back  to 

the  actual  system  xyz  by  (5-  2),  the  relative  velocity  wc 
of  the  core  with  respect  to  the  shell  is 


a  a 

cx2z  '  Fx3y 

b  b 

IX3X  '  cxl2 

_Fxly  "  fx2x 


(5  -  5) 


For  a  rigid  body  the  velocity  w  with  respect  to  inertial 
space  is  obtained  by  (cf.  TNP  I,  p.8): 


V  -t  „ 

|  U2Z  w3y  j 

i  i 

w=wxx=  j;o3x-u)1zj  .  (5-6) 

!  3  i 

[_  uty  "  ^2*  j 

Since  the  mantle  is  rigid,  the  velocity  of  any  particle 
of  the  mantle  with  respect  to  inertial  space  is  given  by  (5-6) 


v  .  =  w 

-man  tie  - 


x  x 


(5-7) 
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For  a  particle  of  the  liquid  core, 
must  be  added  to  (5-6): 

v  =  w  +  w 
—core  —  — c 

Now  the  kinetic  energy  T 


the  relative  velocity  (5-5) 


(5  -  8) 

can  be  evaluated.  We  have 


T 


1  (i 

ea  r  th 

1 

7 


i  i  ( 

vzdM  =  £  j  v  • 


earth 


r  < 


v  »  vdM  + 


V  '  v 


1 


mantle 


( 

v 

core 


T  ,  +  T 

mantle  core 


VdM 

*  VdM 


(5  -  9) 


The  substitution  of  w  »  v^^^  ,  eq.(5  -  7),  gives 


T  »  i  f  |  (  w  •  wdM 

mantle  2  j  j J  —  — 

mantle 


(5  -  10) 


whereas  with  v  =  v  ,  eq.(5-8),  we  get 

-  —core 


1  '  f 

!  -ir  !  i  ( w •  w  +  2 w •  w  +  w»w)dM 

core  c  Jjj  —  —  -  -c  —  c  —  c 

co  re 


(5  -  11) 


The  sum  of  these  two  equations  finally  gives 


1  '  f  < 


l'"  ( 


T  =  t  w-w  dM  +  7T.  (2w*w  +  w  *w  )dM 

2  ;  j  J  -  -  2  j  |  ;  —  — C  — c  — c  ' 

earth  core 


(5  -  12) 
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The  first  term  on  the  right-hand  side  of  (5-12)  is  the 
kinetic  energy  of  a  totally  rigid  earth  since  w  is  given 
by  (5  -  6) .  Thus  ,  by  (4  -  16) , 

jjj  w  •  wdM  =  +  Bu>*  +  (5-13) 

earth 


where 


A  =  jjj  (y2+z2)dM  • 

earth  (5-14) 

B,  C  by  cyclic  permutation  , 

are  the  earth's  principal  moments  of  inertia  as  usual. 

Note  now  that  the  third  term  in  (5-12)  differs  from 
(5  -  13)  by  w  being  replaced  by  wc  .  The  comparison  of 
(5-5)  and  (5-6)  shows  that  this  difference  consists  in  u>. 
being  replaced  by  Xj  ,  and  also  in  the  factors  a/b,  a/c, 
b/c,  etc.  Thus  it  is  not  difficult  to  see  that  the  analogue 
of  (5-13)  is 


'  (  !  2 

w  •  w  dM  =  A  * .  + 
l J )  -c  -c  cA 1 

co  re 


x?  - 


c  *  2 


2 

:X3 


(5  -  15) 


whe  re 


y2dM  + 

.1  )  i  c  '  ; 

core  core 


z 2  d  M 


i  b 


16) 
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3c  and  Cc  following  by  cyclic  permutation  of  x,  y,  z 
and  a  ,  b ,  c 

Similarly  the  second  term  in  (5-12)  becomes 


j  +  Gu>2x2  +  ^3X3  (5_17) 

co  re 


whe  re 


core 


G,  H  by  cyclic 


b 

c 


z  dM 


core 


permutati on . 


(5  -  18) 


The  core  being  a  homogeneous  ellipsoid  of  axes  a,  b,  c  , 
the  integrals  are  easily  evaluated  (actually,  it  suffices  to 
evaluate  one  of  them;  the  others  follow  by  cyclic  permutation), 
and  we  get 


c=|Mc(b2+c2)  , 

F  =  be  . 

5  c 

c  =  k<c2+a2>  • 

G  -  -r-M  c  a  > 
b  c 

(5  -  19) 

c  •  K<a2+b2>  - 

H  =  |m  ab  • 

5  c 

Thus  the  kinetic  energy  (5  -  12)  becomes  finally 
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T  =  +  Bw2  +  Cu)3  +  Acxi  +  Bcx2  +  Ccx3^  + 

+  Fu 1 x 1  +  Gw2x2  +  ^“3X3  •  (5-20) 

Now  we  are  ready  to  apply  Poincar§'s  equations  (4-  13). 
We  obviously  have  6  degrees  of  freedom  (n=6)  .  The  six 

possible  infinitesimal  transformations  are: 

1.  a  rotation  of  the  whole  earth  with  respect  to  iner¬ 
tial  space,  with  angular  velocity  components  Wj,  u>2,  and 
generators  denoted  by  X 1 ,  X2>  X^  and 

2.  a  rotation  of  the  unit  sphere  (5-3)  representing  the 
motion  of  the  core  with  respect  to  the  mantle,  with  angular 
velocity  components  X]*  X2*  X3  and  generators  denoted 


Thus  the  whole  group  relevant  for  the  present  problem 
consists  of  two  independent  rotation  groups.  The  structure 
equations  (3  -  14)  for  rotation  groups,  for  generators  X. 


instead  of 

matri ces  E . 

—  1 

,  give 

[X,. 

li 

X 

tv  V 

II 

1 

-< 

[x2, 

V 

*  X1  • 

tv  V 

-  V  • 

(5  -  21) 

V 

*  X2  , 

tv  V 

11 

1 

-< 

K) 

Note  the  difference  in  sign  due  to  the  fact  that  the  second 

rotation  is  with  respect  to  the  mantle  whereas  the  first  is 

a  rotation  of  the  mantle  with  respect  to  inertial  space.  Any 

rotation  X.  commutes  with  any  rotation  Y.  since  the  two 
1  J 

rotation  groups  are  independent  of  each  other,  whence 
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[  v  v  =  0 


(1  and  j  -  1,  2,  3).  (5-22) 


The  six  quantities  x^  »  x2  »  X3  may  be  identified 

with  Uj,  u>2,  ...  ,  u6  ,  and  X]t  X2  ,  X3 ,  Y,,  Y2,  Y  with 
X,  ,  X  ,  ...»  X6  according  to  sec.  3.  The  corresponding 

constants  of  structure  c...  are  all  0,  1,  or  -1,  bv 
(5-21)  and  (5-22). 


gi  ves 
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torque  L  acts  on  the  whole  earth;  there  is  no  external 
torque  which  would  effect  a  relative  motion  of  the  core  with 
respect  to  the  mantle.  This  relative  motion  is  caused  purely 
by  the  rotation  of  the  mantle  which,  through  the  elliptical 
core-mantle  interface,  acts  on  the  core  through  "inertial 
coupl i ng" . 

Thus  we  have  derived  eqs.  (12-8)  and  (12-9)  of  TNP  I, 
where  their  solution  and  further  implications  are  discussed. 

This  application  already  indicates  the  power  and  useful¬ 
ness  of  Poincarfe's  principle  for  the  mathematical  study  of  the 
earth's  rotation. 


6. 


Elastic  Mantle  and  Liquid  Core 


Kinetic  Ener 


M. 


After  this  digression  to  a  rigid  mantle  we  shall  now 
continue  the  discussion  of  the  solution  of  (Sasao  et  al  . , 
1980)  from  the  point  of  view  of  Poincar§'s  variational  prin¬ 
ciple.  Thus  the  mantle  is  again  considered  elastic. 

The  velocity  vector  v^  of  any  point  of  the  mantle 
with  respect  to  inertial  space  may  be  written  in  the  follow¬ 
ing  way: 

v  =  v  .  =ujxx  +  v=  w  +  v  .  (6-1) 

—  —man  tie  —  —  — m  —  — m  v  ' 


Here 


w  =  w  x  x  (6-2) 

represents  a  rigid  rotation  of  the  coordinate  system  x^, 
x2>  x^  ,  expressed  as  the  vector  product  of  the  angular 
velocity  vector  u  =  (uj  ,u2,Uj)  and  the  position  vector 
*  =  (Xj.x^.x^)  ,  referred  to  body-fixed  axes;  and  _vm 

represents  a  small  residual  velocity  describing  the  rate  of 
elastic  deformation. 

The  meaning  of  the  decomposition  (6-1)  is  easily  under 
stood.  If  the  earth  is  not  rigid,  then  there  is  no  coordinate 
system  at  which  all  particles  of  the  earth  are  at  rest.  Thus 
the  particles  of  the  mantle  move  with  respect  to  our  system 
:< ,  ,  x_ ,  x,  with  velocity  v  which  is  considered  small 
since  it  is  zero  for  a  rigid  body. 

In  the  core  we  similarly  have 


Now  there  is  a  small  rotation 

^  slxX  (6-4) 

of  the  core  with  respect  to  the  mantle,  described  by  the  an¬ 
gular  velocity  vector  x  *  (xi’x2’x3  )  »  superimposed  on 

w  ;  in  addition,  there  are  very  small  residual  motions  vc 
which  describe  deviations  of  the  fluid  core  motion  from  a 
rotation  of  the  core  as  a  whole,  including  effects  of  non¬ 
sphericity  . 

The  decomposition  (6-3)  follows  (Sasao  et  al  .  ,  1980, 

p.167).  It  assumes  that  the  hydrodynamic  core  motion  is  essen 

tially  a  rotation  w  ,  with  small  perturbati ons  v  .  This 

— c 

assumption  is  valid  not  only  for  the  Poincare  model  with  a 
homogeneous  liquid  core  (Poincare,  1910,  sec. 1. 2)  but  very 
generally:  it  corresponds  to  a  toroidal  oscillation  t| 
proportional  to  the  distance  r  from  the  earth's  center, 
which  holds  for  practically  all  models,  even  with  a  hetero¬ 
geneous  core;  cf.  TNP  II,  pp. 135-136. 

The  present  decomposition  (6-3)  is  similar  to  (5-8), 
but  there  are  two  differences.  First,  there  is  an  additional 
term  v  ,  and  second,  wc  in  (6-4)  is  an  exact  rotation, 
whereas  (5-5)  deviates  from  an  exact  rotation  by  terms  on 
the  order  of  the  core  flattening.  These  terms,  together  with 
aeviations  from  the  simple  "Poincar&  motion",  are  now  incor¬ 
porated  into  the  residual  velocity  v 


Hence  (6  -  5)  becomes 


T  = 


1  (fi 


w.wdM  +  7  (2w*w  +  w  »w  ) dM  +  t  +  t  .  (6-1 

—  —  2  —  —  c  — c  — c  m  c 


earth 


core 


(6-2) 


For  the  term  t  we  obtain  from  (6-7)  by  means  of 


t  a  if;  W«V  dM 
m  I - m 


mantle 

i  ( i 


J  ) 

mantle 


(  oi )  •  v  dM 
—  — m 


i  i 


•  jjj  2- (£«*„,)« 

mantle 


( 


-)) 

mantle 


r 


xxv  dM 
—  — m 


(6  -  1 


using  a  well-known  vector  identity  and  the  fact  that 
a  constant  with  respect  to  integration. 
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xxv  dM  =  h 
—  — m  — m 


(6  -  13) 


mantle 


is  nothing  else  than  the  relative  moment  of  inertia  for  the 
mantle;  cf.  TNP  I,  pp. 13-14;  TNP  II,  pp. 110-111.  Substituting 
(6-13)  into  (6-12)  we  have 


t  —  to  -  h  ; 
m  -  — m 


(6-14) 


in  a  similar  way  we  find 


(<£  +  x)  hc 


(6-15) 


where  h_  is  the  relative  angular  momentum  for  the  core. 

It  is  now  of  basic  importance  that  both  hm  and  hc 
can  be  made  zero.  For  the  mantle,  hm  vanishes  if  we  take 
Tisserand  axes  for  the  mantle  as  body  axes  x1x2x^  ;cf. 

TNP  II,  p . 111.  For  the  core,  hc  can  be  made  zero  by  an 
appropriate  definition  of  the  core  rotation  _x  in  (8-3): 
if  h  is  not  zero  for  a  certain  ^  ,  replace  x  by 

x_  +  5^  and  determine  the  three  components  of  6x  by  the  three 


conditions  h_  =  0 

Thus  t  and  t 
m 

that  (6-  11)  reduces  to 


Then  also 


will  vanish. 


can  always  be  considered  zero,  so 


T  =  x-  •  w  •  wdM  +  i  j  f  [  (  2  w  •  w  +  w  _ •  w  )  d H 
2  )  :■ ; —  2  1 J  j  — c  -c  -c 

earth  core 


(6  -  15  ) 


56 


By  manipulations  familiar  from  sec. 5  we  obtain 


T  =  wC.  .u.id .  +  Cf.u.v.  +  icf.x- 


•  J  I  J 


.  .u  .  Y  . 

I J  I AJ 


i  j  x  i  *  j 


(6-  17) 


C.^.  denoting  the  inertia  tensor  for  the  whole  earth,  and 

Cc.  ,  for  the  core  only. 

>  J 

In  view  of  F  =  A  and  similar  relations  following  from 

(5-19),  this  reduces  to  (5-20)  for  a  rigid  mantle,  in  spite 

of  a  slightly  different  definition  of  ^  .  This  provides  a 

check  for  the  decomposition  (6  -  3)  which  treats  wc  as  an 

exact  rotation  (6-4)  and  incorporates  residual  ellipsoidal 

effects  into  v 
— c 

Following  (Munk  and  Macdonald,  1960,  p.37;  TNP  I, 
pp. 15-16)  we  shall  split  up  the  inertia  tensor  as  follows: 


C  =  C  +  c 
—  — o  — 


(6  -  18) 


where 
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L° 

0 

o 

r 

o 

C2  3 

- J 

r*> 

U 

Thus  , 

C 

—  Q 

corresponds 

to  the  model  of  an 

unde  fo  rmed 

earth 

whose 

principal 

axes 

of  inertia  coincide 

wi  th 

the  coordinate 

axes,  which  has  rotational  symmetry  (3=  A)  and  whose  prin¬ 
cipal  moments  of  inertia  A  and  C  are  constant  in  time. 
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The  tensor  £  takes  into  account  the  deviation  of  the  ac¬ 
tual  earth  from  this  simplified  model. 

We  thus  split  up  the  inertia  tensor  C.^.  and  do  the 

same  for  C?.  .  Of  the  components  of  the  residual  inertia 

'  J  c  c 

tensors  we  only  retain  c^,  c23*  ci3  anc*  c23  are 

related  to  nutation  and  polar  motion,  cf.  TNP  I,  sec. 3.  Other 

terms  do  not  influence  these  phenomena  and  can  be  disregarded 

without  harm.  Then  (6-  17)  takes  the  final  form 

T  =  ^(oj^+o)^  )  +  +  C]3W]U3  +  C23a)2u3+  "  20) 

+  +  X*  +  2u2x2  +  X2)  +  7Cc^2w3x3  +  x  3 )  + 

+  c13^1x3  +  U13xi  +  X  |  X3 )  +  c23^w2y3  +  “)3X2  +  X2X3^ 


« 


For  an  ellipsoidal  earth  with  a  solid  mantle  (solid 
earth,  Poincarfe  model),  we  have  a  gravitational  potential 
energy  U  which  depends  on  the  orientation  of  the  body 
with  respect  to  the  directions  to  sun  and  moon  and  hence  on 
the  Euler  angles  4,  0,  <j>  or  the  nonholonomic  variables 
drt  j  ,  dir2,  dir  ^  ;  its  deri  vati  ves  wi  th  respect  to  these  vari¬ 

ables  are  the  torque  components 


all 


cf . sec . 3 . 

For  an  elastic  mantle,  we  have  to  add  a  potential  energy 

U  ,  related  to  the  elastic  deformation,  so  that  the  total 
e  a  — -  ■  -  —  -  - 

potential  energy  is 


U 


U  + 
9 


(7  -  2) 


The  potential  energy  Ued  arises  from  two  causes: 
from  the  elastic  forces  which  are  a  reaction  of  the  earth 
(including  its  liquid  core)  to  the  external  forces,  and  from 
a  change  of  the  gravitational  energy  due  to  the  elastic  de¬ 
formation.  We  thus  have 


where  Ue  represents  the  elastic  energy  in  the  narrow  sense 
and  will  be  called  the  deformation  energy.  Thus  the 

total  potential  energy  becomes 

U  «  U  +  U  ♦  U  .  .  (7-4) 

g  e  d  v  ' 

Because  of  the  smallness  of  U  and  U.  and  because 

e  d 

of  the  smallness  of  the  earth's  flattening  f  ,  the  effect 
of  f  on  Ued  is  of  second  order  and  will  be  neglected. 
Thus,  for  the  purpose  of  computing  U  ,  we  may  use  the 
spherical  approximation,  formally  replacing  the  earth  by  a 
sphere. The  dynamical  effects  of  the  earth's  rotation  on  Ugd 
are  also  of  the  same  second  order  (this  is  easily  understood 
because  the  flattening  is  an  effect  or  the  earth's  rotation) 
and  will  be  neglected.  Thus  again  for  the  purpose  of  compu- 
ting  Ugd  ,  we  formally  consider  the  earth  nonrotating  and 
without  other  kinetical  effects  (with  one  exception  discussed 
below);  this  is  the  static  approximation.  Both  approximations 
are  standard,  from  (Jeffreys,  1949)  to  (Sasao  et  al .  ,  1980); 
only  (Wahr,  1981,  a,  b,  c  >  1982)  is  more  accurate. 

The  gravitational  energy  of  a  material  particle  of  mass 
m  in  a  field  of  potential  V  is  -m.V.  and  that  of  a 

i  i 

system  of  particles  thus 

-Em . V .  .  (7-5) 

The  minus  sign  comes  from  the  tact  that  the  potential  V  in 
geodesy  is  defined  with  opposite  sign  as  compared  to  the 
usage  in  physics.  For  a  continuous  mass  distribution,  (7-5) 
must  be  replaced  by 


dM  denoting  the  element  of  mass. 

Thus  Ud  ,  the  change  of  gravitational  energy  because 
of  deformati on,  becomes 


V  dM 
e 


JJjV 


(7  -  7) 


S 1  denoting  the  deformed  earth  surface.  It  is  essential 
in  this  context  that  the  undeformed  earth  surface  is  consi¬ 
dered  a  sphere  S  ,  so  that  the  deviation  of  S.  from 

o  * 

the  sphere  Sq  represents  the  elastic  deformation  only: 
the  ellipticity  of  the  earth  (the  "equatorial  bulge"  causing 
the  torque  L.  )  has  already  been  incorporatea  into  U  ; 
cf.(7  -  I). 

W  i  th 


U  =  < ( xzcosa  t  +  yzsinct)  (7-8) 

d 

as  usual  (cf.(l-lO)),  eq . ( 7  -  7  )  becomes 


U 


d 


k  c  0  S  o  t 


f  ( 


jj 

s 


xzdM 


<  s i n  •  t  j  j  . yzdM 

s', 


(7-9; 


tne  second  intepral  in  (7-7)  is  zero  because  of  symmetry;  c£ 
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By  the  usual  definition  of  the  products  of  inertia  we  have 


c 


13 


=  -D 

xz 


' 

f 

,  J 

. 

JJJ 


xzdM 


(7  -  10) 


and  similarly  for  ;  cf.  the  analogous  expressions  for 

thecore,  (2  -  12)  and  (2  -  15).  Thus  we  simply  have 


Ud  =  c  1  ^kcoso t  +  c  Ksinot  .  (7-11) 

Using  (1-11)  this  may  be  written 


Ud  =  q(c13W1  +  C23W2}  •  (7  ‘  12) 

This  is  the  deformation  energy. 

The  elastic  energy  Ug  can  be  readily  found  by  the 

following  consideration.  According  to  the  static  approximation, 

external  forces  are  counteracted  by  elastic  forces  in  such  a 

way  that  both  systems  of  forces  are  in  equilibrium.  By  the 

very  basic  principles  of  elasticity,  the  equilibrium  is  reached 

when  the  body  is  deformed  in  such  a  way  that  the  potential 

energy  is  a  minimum.  According  to  the  spherical  approximation, 

the  total  potential  energy  reduces  to  U  ,  (since  U  =0 

ed  g 

for  a  sphere:  if  the  earth  were  a  sphere,  there  would  be  no 
torques  l.  and  no  precession  and  nutation!).  Hence,  for 
equilibrium. 


U  J  =  m i n i m u m  (7-13) 

e  d  v  ' 


j 
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or 


dU  »  0  .  (7-14) 

ed 

By  ( 7  -  3 )  thi s  means 

dU  =  -dl!  .  .  (7-15) 

e  a 

As  far  as  the  external  gravitational  potential  is  con¬ 
cerned,  UJ  is  giyen  by  (7-12).  To  this  we  have  to  add, 
d 

however,  a  term 

+  C23*2}  »  <7  ■  16> 

which  is  to  be  explained  as  follows.  In  the  static  approxi¬ 
mation  used  so  far  we  have  neglected  the  earth's  rotation  as 
expressed  by  complex  number  u  ;  c f . ( 1  -  3 ) .  It  is  known, 
however,  that  the  core  rotation,  as  expressed  by  the  complex 
number  v  ,  can  be  considerably  larger  than  u  .  (In  fact, 
the  table  on  p .  1 2 2  of  TNP  II  shows  that  the  ratio  |v/u!  =  3/e 
can  reach  values  of  200  and  more).  Thus  even  when  we  neglect 
u  ,  we  cannot  neglect  v 

Now  the  incremental  centrifugal  potential  of  the  earth' 
rotation  is,  by  TNP  II,  p.  115, 

2 

-  -Q  e(xzcos,jt  +  yzsinat) 

-  -£(-j  (  xz  +  yz  )  ,  (  7  -  17  ) 

by  (2-3)  and  (2-4).  For  core  rotation,  we  must  replace 


(7  -  18) 


<p  =  -fi(Xlxz  +  y.yz) 

core  XA1  *2  ' 

This  expression  is  now  substituted  for  V  in  (7-6),  the 
integral  being  extended  over  the  core.  The  result  is 

-n(cf3Xi  +  c23x2)  .  (7-19) 

as  the  comparison  of  (7-8)  and  (7-18)  shows:  there  corres¬ 
pond  -nxj  to  <cosat  and  ~n\2  to  <sinat  ,  and  we 
get  cj3  and  for  the  core  since  the  integral  is  exten¬ 

ded  over  the  core. 

Now  we  finally  have  to  add  (7  -  19)  to  (7  -  12)  to  get 
the  expression  for  U  for  the  present  purpose,  denoted 
by  Ud': 

V  '  a(cl3Wl  *  C23W2  '  C*3X1  *  c23*2>  •  (7-20) 

It  may  be  asked  why  U,  ,  and  not  U‘  ,  is  used  in  the 

d  d 

potential  energy  (7-4).  The  reason  is  thav  the  term  (7  -  19) 
is  alredy  incorporated  in  the  kinetic  energy,  so  that  the 
Lagrangian  function  £  =  T  -  U  does  contain  it,  and  it  can 
contain  it  only  once.  (The  reader  may  find  a  similar  reason¬ 
ing  in  (Jeffreys,  1949).) 

For  the  present  purpose,  the  computation  of  the  elastic 
energy  Ue  from  the  equilibrium  condition  (7-13),  does 
however  require  the  use  of  (7-20),  so  that  (7-15)  is  tc  be 
replaced  by 


dU  =  -dll  *  .  (7-21) 

e  a 

By  (7-20)  we  have 

*dV  =  2'-w,dc,  3  *  W2dC23  *  xldC13  +  x2dC23>-  <7  '  22> 

since  the  variations  dUp  and  dU  '  refer  to  variations 

e  d 

of  the  shape  of  the  body,  which  is  expressed  by  changes  dc_k 

and  dcc  in  its  inertia  tensor, 
i  k 

Now  we  shall  use  (2  -  10)  and  (2  -  48)  with  u  =  0 


(static  case): 


c  -  -Dnw  *  0)2V  . 


c  =  -0  w  +  D  v 
c  2  1  22 


(7  -  23) 


Denoting  by  E..  the  elements  of  the  matrix  (D.^)  »  we  have 


a  Enc  +  Ei2cc  • 
v  =  E2,C  *  E22Cc  ■ 


With  w  =  wf  +  iw2  ,  v  =  Xl+  iX2  this  gives 


(7  -  24) 


c 

+  E  cc 

1 1 

13 

12  13 

c 

+  E  cc„ 

2  1 

13 

22  1  3 

and  similarly  for  w 2>  .»2>  c23,  c23 


(7  -  25) 


.  We  substitute  this  into 
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(7-22)  and  take  (7-21)  into  account.  The  result  is 

dUe  =  Q f  1  1  ^ C1  3dC1  3  +  C23dc23  ^ 

+  E12^CUdC13  +  C23dC23} 

+  E21(c13dc13  +  c23dc23) 

+  E22^C1 3dC 1 3  +  C 2 3 d C 2 3 ^  _ 


(7  -  26) 


This  expression  is  a  complete  differential  if  and  only 


E  2 1  =  E  1 2  • 


(7  -  27) 


This,  however,  is  the  case  since  from  the  symmetry  of  the  ma¬ 
trix  D  ,  expressed  by  (2-50),  there  follows  the  symmetry 
of  the  inverse  matrix  E  and  hence  (7-27). 

Thus  (7-26)  represents  the  complete  differential  of  the 
function 


7r<  !  E  11  (  c  13  +  C  23  ^  + 


+  2  E12(c,  3ci 3  +  c  2  3 c  2  3  ^  + 

+  E 2 2 ^ C 1 3  +  C23  )  1  ’ 


(7  -  28) 


which  constitutes  the  desired  elastic  potential. 

Now  the  tidal  potential  energy  related  to  elastic  defor¬ 
mation  is  given  by  (7-3)  as  the  sum  of  (7-28)  and  (7-12): 
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Ued  ■  7nLEl.(CU2+  C232)  +212<C13CU  +  C23C23>  + 

+  E22*c?32  *  c  2  3 2 )]]  +  n  f c  l  3W  l  +  C23W2^  •  <7"29) 

The  great  simplicity  of  this  expression  as  compared  to  analo¬ 
gous  formulas  in  (Jeffreys,  1949)  and  (Jeffreys  and  Vicente, 

1957  )  is  another  indication  of  the  extremely  fortunate  choice 
of  c^,  c2j,  c^.and  c^  as  variables  character!' zi ng  the 
elastic  deformation  of  the  earth;  this  is  the  basic  discovery 
of  (Sasao  et  al . ,  1980)  . 

In  (Moritz,  1982  a)  we  have  found  (7-29)  by  an  ad-hoc 
reasoning  (choosing  Ugd  so  that  the  right  result  comes  out); 
the  present  derivation  thus  provides  a  physical  verification. 

An  independent  check  is  obtained  by  using  a  theorem  con¬ 
cerning  the  potential  energy  of  deformation  given  on  p.173  of 
(Love,  1927):  "The  potential  energy  of  deformation  of  a  body, 
which  is  in  equilibrium  under  a  given  load,  is  equal  to  half 
the  work  done  by  the  external  forces,  acting  through  the  dis¬ 
placements  from  the  unstressed  state  to  the  state  of  equilibrium. 
It  is  not  difficult  to  see  that  this  implies  that  Ug  is  equal 
to  half  of  - U d 1  .  Now  in  fact,  inserting  (7-25)  into  (7-20) 

and  dividing  by  -2  we  get  (7-28). 


I 
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Application  of  Poincarfe's  Principle 


Just  as  in  sec. 5,  we  have  two  rotation  groups,  with 
generators  X ^ ,  X2,  and  Y2,  Y^  •  We  have  used  the 

same  symbols  Xj»  x2»  X3  *•  Y,.  Y2,  Y ^ ,  for  the  core  rotation 
as  in  sec. 5  although  the  conceptual  meaning  is  slightly  differ¬ 
ent:  now  the  core  rotates,  whereas  in  sec. 5  we  had  a  rotation 
of  the  auxiliary  sphere  (5-3). 

In  addition  to  w .  and  x.  ,  the  kinetic  energy  (6-20) 
also  contains  the  variable  products  of  inertia  c]3>  c23’  cn’ 
c\  ;  also  the  potential  energy  (7-29)  depends  on  these  va¬ 
riables.  These  four  quantities  constitute  four  additional 
variables  (four  additional  degrees  of  freedom),  which  describe 
the  elastic  deformation.  They  are  ordinary  (holonomic)  vari¬ 
ables  q  ,  qg,  q^,  q]Q  ,  so  that  (4-3)  holds  for  them. 

(This  also  fits  into  the  group  theoretic  scheme,  with  * 
instead  of  q.  for  7  s  i  s  10  ,  the  correspondi ng  subgroup 

beinq  Abelian  with  zero  c..,  .) 

1  jk 

The  Poincare  equations  (5-23)  and  (5-24)  remain  finally 
the  same  since  we  have  two  independent  rotation  groups  as  in 
sec. 5.  In  addition  to  these  six  equations  we  have  the  equations 
3  £ / 3  q .  =  0  for  i  =  7,8,9,10  ,  which  give 


The  10  equations  (5  -  23),  (5  -  24),  (8  -  1)  and  (8-2) 
relate  and  determine  the  10  quantities  u  1  ,ui2  ,x  j  »x2  ‘Xj  »c  1  3  * 

c23  ’C1 3  *  ^2  3 

We  note  that  the  torque  components  L.  are  the  same 
as  in  (5-23),  namely  purely  gravitational;  in  fact,  by  (7-1) 
and  ( 7  -  2 )  we  have 


since  (7-29)  does  not  depend  on  rotation  so  that  a Ue d / ;  =  0. 

On  the  other  hand,  the  right-hand  sides  of  (8-1)  and  (8-2)depend 
only  on  the  energy  of  elastic  deformation  since 

all  _  3Ued 
3C1  3  '  3 c  1  3  ’ 

as  U  does  not  contain  c,,  ,  and  similarly  for  the  other 

g  1  j 

elastic  variables. 

Let  us  now  substitute  (6-20)  into  the  third  equation 
of  (  5  -  24 ) .  This  gives 
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Now  we  have 

u1 »  u2  <<  Wj  ;  C1 3 ’  c23  <<  A,  C  ,  (8-4) 

(if  we  take  absolute  values),  and  similarly  for  the  core; 
because  of  the  smallness  of  the  core  rotation,  also  y 3 
will  be  smal 1 .  Thus 

ur  10  2  ’  xr  x2  ’  x3  ’  C 1  3  ’  °2  3  ’  CU’  C23  (8~5) 

can  be  regarded  as  small  quantities  of  the  first  order;  their 
squares  and  products  are  then  of  the  second  order  and  will  be 
negl ected . 

Then,  and  since  Cc  is  constant  by  definition,  (8-3) 
reduces  to 

K<Cc"3  *  Cc*3>  ■  °  '  (8-6) 

If  we  treat  the  third  equation  of  (5-23)  the  same  way,  sub¬ 
stituting  (6-20)  and  neglecting  second-order  quantities,  we 
ge  t 


cft(  Ctu3  +  Ccx3}  =  L3  =  0 


(8  -  7) 


since  =  0  because  of  rotational  symmetry.  The  subtraction 

of  (8-6)  from  (8-7)  gives,  C  and  Cc  being  constant  by 
def ini ti on , 


{jw 

(C  '  Cc 


0 
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so  that  uij  =  const.  We  may  put 


w  3  =  ft  »  (8*8) 

identifying  <j3  with  the  average  rotational  velocity  ft  of 
the  earth;  this  is  possible  since  u  and  u2  are  very  small. 

Now  (8-6)  yields  dX3/dt  =  0  ,  so  that  =  const. 

The  simplest  choice  for  this  constant  is  zero,  giving 


x  3  =  0  (8-9) 

This  is  the  solution  we  shall  take. 

If  we  subs ti tute  ( 6  -  20 )  into  the  first  two  equations 
of  (5-23),  neglect  second-order  quantities  and  use  (8-8) 
and  ( 8  -  9 ) ,  we  obtai n 


Ai,  +  (C  -  A)flu2  +  a c(x,  -  nXz)  +  n(c13  -  nc23)  =  L1  , 

(8-10) 

A(L2  -  (C-A)nUl  +  Ac(x2  +  nXl)  +  ii(c23  +  nc13)  =  l2  . 

If  we  do  the  same  with  the  first  two  equations  of  (5-24), 
we  get 


V,  +  Ac*l  "  Vx2  +  =  °  ’ 

Ac"2  +  Ac  x  2  +  Vx,  +  8^5  "  0  ’ 


(8-11) 


Let  us  finally  substitute  (6-20)  and  (7-29)  into 
(8-1)  and  (8-2),  again  using  ^  =  ft  and  x3  =  0  .  The 

result  is 


wi  a  Enci 3  +  E 1  2 C  1  3  +  W1  ’ 


w 2  *  E  1  1  C 2  3  +  E 1  2 C 2  3  +  w2  ’ 


(8  -  12) 


*1  =  E  1  2 C  1  3  +  E22CU  ’ 
x2  “  E  1  2C 2 3  +  E 22° 2 3  ' 


(8  -  13) 


Now  the  complex  combination  of  (8-10)  and  (8-11)  gives 
Au  -  i ( C  -  A)qu  +  A^ ( v  +  i av  )  +  a(c+inc)  =  L  ,  (8-14) 


Aii  +  Av+iCr.v+nC  =  0  , 

c  c  c  c 


(8  -  15) 


and  the  complex  combination  of  (8-12)  and  (8-13)  becomes 


u  =  Ei;c  *  Eizcc  *  w  ' 


»  «  E|2c  ♦  E22cc  . 


(8  -  16) 


3y  means  of  the  matrix  r D . . ]  inverse  to  the  matrix  [E 

'  J  i  j 

the  last  two  equations  may  be  written 


c  -  0 1 , ( u-w)  +  0 1 ^ v 


(8  -  17) 


cc  =  D , 2 ( u-w )  +  D22v 


(8  -  18) 


4 
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Our  equations  (8-14),  (8-15),  (8-17),  and  (8-18) 
are  identical  to  the  basic  SOS  equations  (1-1),  (1-2),  (1-6), 
and  ( 1  -  7 )  of  sec . 1 . 

Concluding  remarks.The  present  variational  method  fur¬ 
nishes  probably  the  most  direct  derivation  of  equations  (1~1) 
and  (1-2),  or  (8-14)  and  (8-15).  This  is  particularly  true 
for  (1-2)  which  has  been  derived  by  Sasao  et  al.(  1980)  using 
both  the  decomposition  (6-3)  and  the  hydrodynani c  equations. 

The  present  approach  uses  only  (6-3)  and  provides  a  unified 
deduction  of  (1-1)  and  (1-2)  which  also  explains  the  similar¬ 
ity  of  these  two  equations  that  in  the  complicated  hydrodyna- 
mical  derivation  comes  out  almost  as  a  miracle.  In  fact,  the 
similarity  of  (1-1)  and  (1-2)  is  now  seen  to  be  due  to  the 
fact  that  each  equation  essentially  reflects  the  action  of  a 
rotation  group. 

On  the  other  hand,  equations  (1-6)  and  (1-7)  are  better 
obtained  directly,  in  the  way  described  in  sec. 2.  In  fact,  the 
derivation  of  the  elastic  energy  given  in  sec. 7  already  uses 
(1-6)  and  ( 1  -  7)  . 

To  summarize:  of  the  four  SOS  equations,  (1-1)  and 
(1-2)  are  obtained  more  easily  from  the  variational  method, 
whereas  (1-6)  and  (1-7)  are  better  derived  directly.  There¬ 
fore,  in  the  first  presentation  of  the  method  described  here 
(Moritz,  1982a)  we  have  followed  a  "hybrid"  approach,  consi¬ 
dering  (1-6)  and  (1-7)  given  (from  a  direct  derivation), 
and  choosing  the  elastic  deformation  energy  Ugd  ,  eq.(7  -  29), 
simply  in  such  a  way  that  the  variational  principle  not  only 
furnishes  the  new  equations  (1-1)  and  (1-2),  but  also  correct¬ 
ly  reproduces  the  given  equations  (1-6)  and  (1-7) 

This  ad-hoc  choice  of  U  was  frankly  pragmatic  and 
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motivated  by  reasons  of  mathematical  simplicity.  It  has  been 
criticized  for  providing  no  physical  interpretation  of  Ugd  . 
Therefore,  we  have,  in  the  present  sec. 7,  derived  Ued  by 
physical  considerations  in  such  a  way  that  the  physical  mean¬ 
ing  becomes  evident.  See,  in  particular,  the  decomposition 
(7-3)  of  the  elastic  deformation  energy  Ued  as  the  sum  of 
an  elastic  energy  Ug  in  a  narrower  sense,  and  a  deformation 
energy  Ud  .  The  latter,  Ud  ,  could  be  determined  immediate¬ 
ly,  and  then  Ug  was  obtained  as  the  elastic  equilibrium 
response  to  Ud  .  Still,  eqs .  (1-6)  and  (1-7)  served  as 

a  basis  also  for  the  physical  derivation  of  U  .  in  sec. 7. 

e  a 

Thus  there  is  no  doubt  that  even  the  use  of  the  present 
variational  principle  does  not  supersede  a  direct  derivation 
of  (1-6)  and  (1-7).  Also,  eq.(l-l)  could  have  been  found 
rather  easily  by  other  considerations,  using  the  Euler  -  Liou- 
ville  equation  (TNP  II,  sec. 10).  There  remains  (1-2),  for 
which  the  present  variational  principle  furnishes  indeed  a 
derivation  of  incomparable  simplicity. 

But  does  this  fact  alone  justify  the  use  of  the  whole 
mathematical  machinery  of  Poincare's  variational  principle? 

The  answer  is  the  same  as  in  many  other  applications  of  ana¬ 
lytical  mechanics:  the  results  could  also  be  found  by  element¬ 
ary  Newtonian  methods,  but  it  is  the  application  of  Laqrangian 
or  Hamiltonian  methods  which  makes  us  better  understand  the 
mathematical  and  physical  structure. 

It  may  be  helpful  to  again  summarize  the  basic  logical 
structure  of  the  present  approach  which  may  have  been  obscured 
by  the  computational  details.  We  start  from  (1-6)  and  (1-7) 
as  deri ved  in  sec  .  2  : 
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c  =  Dn  (u-w)  +  0  2v  , 
cc  *  D12(U'W)  +  °22V  * 


Then  the  kinetic  energy  (6-17) 


ttC  .  .  .  Ul  .  +  C?.UJ.Y.  +  TJ-C  ?  .  X  .  X  • 

2  i  j  i  j  t  j  i  A  j  2  i  j  A  i  A  j 


used  in  the  form  (6-20),  together  with  the  gravitational 
energy  ,  is  substituted  in  Poincare's  equation  (4-14), 


±  fll_l  +  c  31  -  l 

dt  [a.o.j  ci  jk  ”  Li  * 

or  more  explicitly  (5-23)  and  (5-24),  which  directly  gives 
(1-1)  and  (1-2),  or  (8-14)  and  (8-15): 


Au  -  i(C-A)au  +  Ac(v-Hqv)  +  n(c  +  ir,c)  =  L  , 

A  ii  +  A  v  +  iC  qv  +  :iC  =  0  . 
c  c  c  c 

The  gravitational  energy  U  is  not  used  directly,  but 

only  through  the  well-known  expressions  (1-  12)  for  the  torque 

components  L.  .  As  we  have  seen  above  (after  eq.(8-2)  ), 

we  nave  3l)  /air.  =0  ,  so  that  the  elastic  deformation 

e  a  i 

energy  does  not  contribute  to  the  torque.  Hence  the  simple 
expressions  (1-12)  hold  independently  of  the  internal  con¬ 
stitution  of  the  oody. 

We  also  see  that  the  left-hand  sides  of  both  equations 
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(1-1)  and  (1-2)  only  depend  on  the  kinetic  energy  T 
through  Poincar£'s  equation  (4-14),  no  Euler  -  Liouville 
condition  and  no  hydrodynamics  are  needed! 

So  far  everything  has  been  rather  analogous  to  the  simple 
Poincar§  model  (rigid  mantle  and  homogeneous  core)  as  out¬ 
lined  in  sec. 5.  To  take  care  of  elasticity,  we  have  only  to 
add  a  term  Ued  to  the  gravitational  potential  energy  to  get 
the  total  potential  energy 


U  =  U  +  U  , 
g  e  d 


Ugd  is  chosen  in  such  a  way  that  Poincare's  principle, 
through  (8-1)  and  (8-2)  (which  do  not  depend  on  Ug  !) 
gives  equations  (1-6)  and  (1-7);  the  expression  (7-29), 


ed 


Ell^ci3  +C23  )  + 2E12(c13C13+C23C23^  + 


+  E22(c)3  +c23  )  |  +  fi(c13w1+c23w2)  * 


is  determined  uniquely  (up  to  an  additive  constant)  by  this 
condition. 

The  present  variational  approach  may  be  compared  to  that 
of  Jeffreys  (1949)  and  Jeffreys  and  Vincente  (  1957  ).  Common 
to  all  three  approaches  is  the  reduction  of  a  problem  of  con¬ 
tinuum  mechanics  (infinitely  many  degrees  of  freedom)  to  a 
problem  of  analytical  mechanics  with  only  a  finite  number 
(here,  10)  of  degrees  of  freedom.  This  is  brought  about  by 
restricting  ourselves  to  second-degree  harmonic  perturbations 
and  disregarding  all  the  -infinitely  many  other  harmonics. 

The  essentially  greater  simplicity  of  the  present  method 


is  brought  about  by  the  use  of  anholonomic  Poincar§'s  vari¬ 
ables  related  to  the  two  rotation  groups  considered,  and  of 
the  correspondi ng  anqular  velocities  w.  and  x;  (Jeffreys 
uses  holonomic  rotational  variables),  as  well  as  by  the  for¬ 
tunate  choice  of  ^3*  c23’  cn*  an<*  c23  as  e^ast''c  variables 
due  to  Sasao,  Okubo  and  Saito. 

The  resulting  SOS  equations  (1-1),  ( 1  -  2  ) ,  (1-6) ,( 1  -  7  ) 
are  equivalent  to  Molodensky's  equations  both  with  respect  to 
accuracy  (spherical  approximation  for  the  mantle)  and  to 
applicability  to  rather  general  earth  models  (heterogeneous 
mantle,  heterogeneous  core,  even  elastic  inner  core).  They  not 
only  give  the  simplest  formulation  of  Molodensky's  liquid- 
core  problem,  but  also  clearly  show  its  logical  structure. 

Their  practical  usefulness  will  be  seen  in  the  following 
sections. 


PART  C 


EXPRESSIONS  FOR  NUTATION  AND  POLAR  MOTION  INCLUDING 


LIQUID-CORE  EFFECTS 


Transformation  of  Eq .  (1-1 


This  equation 


Au  -  i ( C-A )nu  +  A  ( v  +  i qv  )  +  fi(c+iftc) 


(9  -  1) 


will  be  basic  for  the  investigation  of  the  movement  of  vari¬ 
ous  axes  (rotation  axis,  angular  momentum  axis,  figure  axis, 
etc.)  for  the  Mol odensky-SOS  model  to  be  performed  in  Part  C. 


Let  us  reca 


the  notati ons  : 


denote  the 
I  C  those 


moments  of  inertia  of  the  whole  earth,  A  and  C  those 

c  c 

of  the  core,  a  denotes  the  (average)  angular  velocity  of  the 
earth,  u  and  v  are  the  complex  numbers  describing  small 
rotations  (  u  of  the  whole  earth  with  respect  to  space, 

v  of  the  core  with  respect  to  the  mantle),  and  L  is  a 
complex  number  describing  the  lunisolar  torque;  these  complex 


numbers  are  defined  by  (1-3)  and  (1-5).  The  complex  number 
c  combines  anomalous  elements  of  the  inertia  tensor,  given 
by  (1-4).  It  is  related  to  u  by  (2-8) 

c  =  -jG"1ka5fl(w-u)  ,  (9-2) 

where  G  is  the  gravitational  constant,  a  can  be  considered 
the  radius  of  a  spherical  earth,  k  denotes  the  potential 
Love  number  depending  on  frequency  (we  keep  in  mind  the  possi¬ 
bility  of  the  decomposition  (2-9)  but  shall  not  use  it), 
and  w  is  related  to  the  lunisolar  (tidal)  potential  by  (1-11) 
Thus  ^he  anomalous  elements  of  the  inertia  tensor,  described 
by  c  ,  depend  on  u  (rotational  deformation)  and  w 
(tidal  deformation);  cf.  TNP  I,  p.  27. 

We  seek  solutions  that  are  functions  of  time  t  of  form 

(9  -  3) 

with  real  constants  u  ,  v  ,  taking 

o  o 

w  =  w^e 1 0  c  ;  (9  -  4  ) 

this  is  in  agreement  with  (1-11): 

w  =  a  "  1  <  e  ‘  '3 1  (  9  -  5  ) 


u  =  u  e 
o 


i  a  t 


v  =  v  e 
o 


a  t 


These  particular  solutions  have  the  well-known  advantage  that 
differentiation  is  simply  mul ti pi i ca ti ng  by  io  : 
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u  = 


du 

3t 


=  iju  e 
o 


( a  t 


=  lull 


(9  -  6) 


They  depend  on  the  frequency  a  ;  but,  our  equations  being 
linear,  we  can  obtain  a  rather  general  solution  by  adding  the 
contri buti ons  from  different  frequencies.  On  substituting 
(9-2),  assuming  (9-3)  and  considering  (9-6),  the  basic 
equation  (9-1)  becomes  after  dividing  by  i 


fAc  “  (  C-A  )ft 


u+  A  (a+n)v  -  i-G"  1  ka5n2  ( w-u  )  (  0  +  S2 ;  =  -  i  L 

C  J 

(9  -  7) 

It  is  now  convenient  to  introduce  the  dimensionless 
constant 


k 

s 


3G( C-A ) 
,5  2 


(9  -  8) 


called  secular  Love  number  (TNP  I,  p.  21),  and  express  the 
torque  L  in  terms  of  the  complex  potential  coefficient  w 
by 

L  =  -  i  ( C-A )qw  ;  (9-9) 

this  is  an  immediate  consequence  of  (1-11)  and  (1-12)  (put¬ 
ting  a  =  1  ,  which  is  of  no  relevance).  Thus  (9  -  7)  becomes 


j^Aa  -  (C-A)r,  +  ■£■( a+0  )  (  C-A )J  u  +  Ac(o+n  )  v 
-  -  <C-A)[l  -  • 


(9-10) 


We  now  introduce  the  Euler  frequency  defined  by 


C  —  A 

~ra 


(9  -  11) 


cf.  TNP  I,  p.  10,  and  divide  by  A  .  The  result  is 

A 


a-  a ^  + 


I  (J  rfl  “ 

(o  +  ft)^— Ju  +  ^j~(0+ft) 
> ' 


V  = 


P,  k  o  +  n 
■^l1  ■  IT  — 


(9  -  12) 


In  agreement  with  TNP  I,  p.  121,  we  now  introduce  the 
tidal  frequency  by 

a  =  -ui.  ;  (9-13) 

then 

o.  +  c  =  -U  .-a)  =  -coj  .  (9-14) 

j  j 


whe  re 


is  the  corresponding  nutational  frequency  (TNP  I,  p.  36).  Thus 
we  can  put 


ft  +  a  =  “  Aw  ,  o  =  -(n+Aw) 


(9  -  16) 


omitting  the  subscript.  Dividing  by  -n  we  then  get 


1  +  _i  +  Au  +  k - E  L  +  c  Au 

n  n  k  2  An 

S  46  I 


_JL  i  +  |L_  w 

~  1  r  a  w  • 


(9  -  17) 


Since  the  usual  tidal  frequencies  w .  are  very  close 
to  n  ,  the  angular  motions  of  sun  and  moon  being  so  much 
smaller  than  the  earth's  rotation,  we  have 


(9  -  18) 


Similarly, 


fE  •  1  j 

TU7J  '  ’ 


(9  -  19) 


so  that  the  last  term  between  parentheses  in  the  factor  of 


u  in  (9-  17)  will  be  of  second  order  and  can  be  neglected. 
Thus  we  finally  obtain 


1  + 
l 


l± 

n 


n 


— ]w 

k$  q  j 


A  A 
_£ 

A  n  v 


(9  -  20) 


This  equation  expresses  u  in  terms  of  the  tidal  poten 
t i a  1  characteri zed  by  w  and  of  the  core  rotation  v  .  Let 
us  recall  what  u  means.  It  is  defined  by 

u  =  +  ioj2  (9  -  21 ) 

as  the  complex  combination  of  the  x  and  y  components  of 
the  rotation  vector  ^  .  I f  we  divide  these  components  by 
the  length  a  of  this  vector,  we  get  the  components  of  the 
unit  vector  of  the  (instantaneous)  axis  of  rotation.  The  x 
and  y  components  of  this  vector,  m1  and  m£  ,  can  by  com¬ 
bined  into  a  complex  number 


m  =  mJ  +  im2  (9-22) 

which  describes  the  polar  motion  of  the  rotation  axis  R  : 


PR 


(9-23) 


the  notations  are  the  same  as  in  TNP  I  (cf.  pp.  17  and  93) 
and  in  TNP  II  (cf.  pp.  128  and  138): 
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m,  =  eCOSot 

m_  =  e  s  i  na  t 
2  * 


o  r 


m  =  ee 


i  a  t 


(9  -  24) 


(9  -  25) 


so  that  e  is  the  amplitude  of  polar  motion  for  the  frequency 
a  under  consideration.  Then 


u  =  fie  e 


i  o  t 


We  likewise  put 


v  =  sue 


0 1 


w  =  five 


i  o  t 


(9  -  26) 


(9  -  27) 
(9-28) 


The  numbers  e ,  c  ,  and  t  are  dimensionless  real  constants. 

The  substitution  of  (  9  -  26  ),  (  9  -  27  ),  and  (9-23)  into 
( 9  -  20 )  gi ves 
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0  E  [ ,  kAojl  ^  c  A 

IT, 1  +  IT  “lT  -  X  Xr- 

s  > 


(9-29) 


For  a  rigid  earth  without  liquid  core  (k=0,  A.-  0)  this  gives 
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e 
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-1 


Aw ! 

n  i 


(9-30) 


This  equation  gives  the  amplitude  of  polar  motion  forced  by 

a  tidal  potential  of  frequency  o  and  amplitude  f  .  In 

terms  of  this  c  we  may  write  (9-29)  as 

o 

'  *  ‘o(l  *  TT  if]  -  ‘TTT^  <9-31> 


where  we  have  put 

*  *  (3  *  T  *  Tj"  1  1  '  (9-32> 

In  fact,  in  view  of  (9-18)  and  (9-19),  the  factor  a 
differs  only  very  little  from  unity. 

Special  cases.  For  the  rigid  earth  (k=0,  c=0),  eq.(9-31) 
reduces  to  (9-30),  which  we  can  also  write,  using  (9-15): 

a  E  aE  (9-33) 

£  —  — — . T  =  —  -  J  9 

O  Q  +  A w  +  g  r  u  ,+oP 

E  J  E 

which  is  the  same  as  the  amplitude  of  TNP  I,  eq.(ll-45),  with 

i  =  iB  .  . 

J 

For  the  elastic  earth  without  fluid  core  (k=0.3  ,  c=0), 
eq . ( 9  -  31 ) ,  together  with  (9-33),  gives 


i  j.  k  Aw  °  E 

1  +  F?  77+^7T  • 

s  J  J  E 


(9  -  34) 


On  the  other  hand,  we  have  McClure's  expression  (TNP  I,  p.36), 
apart  from  the  factor  i  (see  next  section): 


<  .B. 
J  J 


t.  Aw  . 

1  +  t)- - »- 

C-A  J  ks  3  g 

"T~  V°c  ,  +  JL  I E  j 

k  q 
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(9  -  35) 


whe  re 


k  j  c  J  E 

I  +  —  -i 

k  sj 


(9  -  36) 


is  the  Chandler  frequency.  Are  (9-34)  and  (9-35)  equivalent? 
We  have  by  (9-11), 


C-A  a 

A  w  .  +  o , 


.  .  +G 

J  C 


so  that  (9-34)  and  (9-35)  will  be  equal  provided 


(9-37) 


+  o , 


fU) 


+  IT  -fi 

c  v 


-J-  -  1 


-  (jj 


Aw  . 


+  ir 

s 


(9-38) 


which  differs  from  iU.+o£  only  by  the  last  term  which  is 

small  of  second  order  since  Aw.  and  a,,  are  both  small,  b' 

J  E 

(9-13)  and  (9-19),  and  which  we  can  neglect. 

For  the  Poincare  model  (rigid  mantle,  liquid  core)  we 
have  k=0  and  c /  0,  so  that  (9-31)  reduces  to 

A  . 

,  c  iU 

■  -o  ■  'X  IT 


(9-39) 


with  \  =  1  .  Since  au/£2  is  very  small  for  tidal  pertur¬ 

bations,  this  is  almost  the  same  as  (9-31). 

From  this  it  would  be  tempting  to  conclude  that  the 
effect  of  elasticity  on  e  is  very  small,  at  least  for  nearly- 
diurnal  frequencies  which  we  have  been  considering  in  this 
section  (this  is  plainly  not  true  for  the  Chandler  period  which 
is  lengthened  due  to  elasticity  from  305  days  for  k=0  to 
403  days  for  k=0 . 3  ).  Thus  we  might  think  that  the  simple 

Poincare  model  (rigid  mantle)  might  be  used  to  calculate 
liquid-core  effects  on  e  ,  still  getting  basically  the  same 
numerical  results  as  for  the  Molodensky  model  with  an  elastic 
mantl e  . 

Such  a  conclusion  would  be  wrong,  however,  as  the  com¬ 
parison  of  the  tables  in  TNP  I,  pp  122-123, and  TNP  II,  p.  122, 
shows.  The  reason  is  that,  though  the  direct  effect  of  k 
on  (9-31)  is  negligible,  there  is  also  an  i ndi rect  effect : 
elasticity  affects  core  rotation,  hence  r  is  different  for 
a  rigid  and  an  elastic  mantle,  and  elasticity  enters  indirectly 
in  ( 9  -  31 )  th  ough  q 

This  shows  that  (9  -  31)  shows  very  well  the  effect  of 
mantle  elasticity  and  core  rotation  if  we  assume  that  the  core 
rotation  parameter  r  is  known.  If  only  the  luni solar  poten¬ 
tial  (coefficient  r  )  is  given,  then  we  must  proceed  differ¬ 
ently,  in  a  way  to  be  outlined  now. 

Solution  for  u  and  v  in  terms  of  w  .  In  addition 
to  (1-1),  let  us  now  also  consider  (1-2): 

Au  +  Av  +  iC;.  v  +  .^c  =  0.  (9-40) 

c  c  c  c  ' 

We  express  c  ^  in  terms  of  u  and  w  by  (2-41)  and  sub- 
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stitute  into  (9-40).  With  (9  -  26),  (9  -  27),  and  (9-28), 
differentiation  is  replaced  by  multiplication  by  i a  ,  accord¬ 
ing  to  (9-6).  This  reduces  the  differential  equation  (9-40) 
into  a  linear  equation  of  the  form 

a2lu  +  a22v  =  C2W  •  ( 9  -  41) 

This  is  in  complete  analogy  to  the  procedure  by  which 
we  have  obtained  (9  -  17)  from  (9  -  1);  we  may  also  write  (9  -  17) 
in  the  form 

anu  +  a]2v  =  Cjw  .  (9  -  42) 

The  solution  of  the  two  linear  equations  (9-42)  and 
(9-41)  for  u  and  v  (say,  by  means  of  determinants)  gives 

u  =  b  w 

v  =  b2w 

Since  it  is  not  difficult  to  see  that  the  coefficients  a.. 

'  j 

and  c.  are  all  real,  also  the  coefficients  b.  will  be  rea  i  . 

i  i 

In  terms  of  the  real  constants  e  ,  <; ,  t  we  have  the  equations 
air:+ai2c=clT 

(9  -  43) 

a  2  1 e  +  a22;  *  V 

with  the  solution 


0  . 


(9-44) 


with  the  same  coefficients  as  above. 


Thus  both  e  and  q  can  be  expressed  in  terms  of  t 
only.  The  computation  is  easy,  but  the  resulting  expressions 
are  not  very  elegant.  Therefore  they  will  not  be  given  here. 

The  corresponding  calculations  for  the  Poincare  model  (k=0) 
can  be  found  in  TN?  I,  pp.  120-122. 

This  reduction  of  a  system  of  differential  equations 
to  a  system  of  linear  equations  (9-43)  has  been  made  possible 
by  using  exponentials  (9-3)  to  (9-5).  This  is  nothing  else 
than  a  trans format i on  from  the  time  domain  to  the  frequency 
domain  which  is  known  to  simplify  matters  on  such  occasions. 

As  a  matter  of  fact,  the  coefficients  a.,  and  c.  ,  and 
hence  b 


a  .  . 

1 J  1 

in  (9-44),  will  depend  on  the  frequency 
Relation  to  the  Molodensky  coefficient  3  .  In  TNP  II, 


C  -  -U) 


p.56,  we  have,  for  the  Poincare  model,  found  the  relation 


v  =  Cl 
o 


(3-e) 


(9  -  45) 


whence  ,  by  ( 9  -  27  )  , 


(9  -  46) 


0  r 


(9  -  47) 


3  being  the  Mo’odensky 
interesting  cases  we  have 
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•>>._  ,  so  that  a  pprox  i  ma  te  I  y 
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It  can  be  shown  that  these  relations  also  hold  for  an  elastic 
mantle.  In  fact,  substituting  (9-47)  into  eq. (10-41)  of  TNP  II 
it  is  not  difficult,  though  somewhat  laborious,  to  derive 
(9-29).  Substituting  (9-26)  and 

v  =  fl(s-£)e'at 


into  the  second  SOS  equation  (9-40)  and  taking  a  homogeneous 
core,  one  obtains  eq.  (9-34)  of  TNP  II.  (Hint:  show  that 


c 

c 


4ttG 

TTpe 


2k5  _ i o t 
b  ye 


» 


notations  as  in  TNP  II,  sec. 2.) 


10.  The  Kinematic  Axes 


In  this  section  we  shall  investigate  the  motion  of  the 
rotation  axis,  the  figure  axis,  and  the  "mean  Liouville  figure 
axis".  These  axes  will  be  called  kinematic  axes,  in  contrast 
to  the  angular  momentum  axis  which,  being  dynamically  defined, 
behaves  differently  and  will  be  considered  in  the  next  section. 

Reference  frame.  As  in  the  previous  reports  TNP  I  and 
TNP  II,  we  shall  use  two  basic  reference  frames,  the  "nutation 
frame"  x°x°x°  which  is  the  natural  frame  for  describing 
nutation,  and  the  "body  frame"  x^x^  to  which  polar  motion 
refers;  cf.  TNP  II,  pp.5-6  and  140-143. 

The  nutation  frame  x°  is  connected  to  the  inertial 
system  in  a  prescribed  way:  the  x°  axis  has  a  fixed  direction 
in  inertial  space,  and  the  system  x°x°x°  rotates  with  con¬ 
stant  angular  velocity  n  around  the  x°  axis.  Being  space- 
fixed,  the  axis  is  a  natural  reference  for  nutation. 

The  body  frame  x.  represents  a  system  of  Liouville 
axes  for  the  mantle  (TNP  II,  pp. 140-143)  with  respect  to  which 
the  mantle,  in  the  absence  of  elastic  deformati ons ,  would  be 
at  rest.  The  x^  =z  axis  represents  the  figure  axis  of  the 
undeformed  ellipsoidal  earth  and  is  the  "mean  Liouville  fi¬ 
gure  axis"  mentioned  above.  It  is  the  natural  origin  of  polar 
motion  and  corresponds  to  the  point  0  in  Fig.  6.1  on  p.38 
of  TNP  I. 

The  basic  trick.  These  two  frames  are  related  by  a  small 
rotation  (TNP  II,  p.  6): 


is  an  "infinitesimal  rotation  matrix".  In  terms  of  the  vector 


J  = 


r°7 


(10-3) 


this  may  also  be  written 


x  =  -  e 


(10-4) 


the  cross  denoting  the  vector  product  as  usual. 

Now  these  relations  between  coordinate  systems  are 
geometrically  the  same  as  the  formulas  for  an i nf i ni tesimal 
rotation  of  a  rigid  body.  In  fact,  the  body  axes  x i x 2 x 3 
can  be  assumed  to  be  fixed  with  respect  to  an  undeformed,  or 
ri gi d, f i cti ti ous  "reference  earth".  Thus  we  may  use,  with 
some  care,  the  simple  formulas  for  the  rotation  of  a  rigid 
body  (TNP  I,  secs. 10  and  11)  even  to  describe  nutation  and 
polar  motion  for  an  elastic  earth  with  or  without  a  liquid 
core.  This  is  the  basic  trick  which  will  be  used  here,  follow 
ing  Wahr(1981c)  and  indications  in  TNP  II,  pp. 136-138. 

Basic  relations.  Denote  the  unit  vector  of  the  rotation 


axis  R  by  e  p  ,  the  unit  vector  of  the  instantaneous  figure 
axis  F  (which  is  defined  as  the  axis  of  maximum  inertia 
for  the  deformed  earth)  by  £  p  ,  and  the  "mean  Liouville 
figure  axis"  (which  is  nothing  else  than  the  x3*z  axis  of  the 
the  body  frame)  by  e?  .  Denote  further  the  unit  vector  of 
the  x°  axis  (of  the  nutation  frame)  by  e  ^  .  Then  polar 
motion  and  nutation  of  the  various  axes  can  by  characterized 
by  the  vectors 
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cf.  Fig.  10.1. 


FIGURE  10.1.  Polar  motion  of  the  rotation  axis 


Since  all  vectors  eR,  ep,  £z,  and  e^  are  almost 
parallel  to  each  other,  the  polar  motion  vectors  £  and  nu¬ 
tation  vectors  £  will  be  small  and  only  have  a  negligible 
component  along  the  x^  axis  (Fig,  10.1).  Thus  we  can  put 


"pr 

£  *  P2 
0 


n  =  f  n. 


(10-6) 


and  characterize  polar  motion  and  nutation  by  the  complex 
numbers 


P  *  P,  +  1 P2 


n  =  n1  +  i n2 


(10-7) 


instead  of  the  three-dimensional  vectors  £  and  £  .  More 
details  about  this  can  be  found  in  TNP  I,  sec. 11. 

First  we  note  the  basic  relation  between  nutation  and 
polar  motion:  from  (10-5)  there  follows  for  any  axis 


p  =  £  +  (e 3~ez )  =  £  -  £ 


or  for  the  correspond! ng  complex  numbers. 


p  =  n  -  n. 


(10-8) 


n  =  p  +  n. 


n  may  be  expressed  by  the  sennd  equation  of  (11-37)  of 
TNP  I,  noting  that  the  figure  axis  F  for  the  rigid  "refer¬ 
ence  earth"  is  the  z-axis: 


e  being  defined  by  the  complex  number 


9  “  9 1  +  1 92 


(10  -  10) 


combining  the  first  two  components  of  the  vector  (10-3).  The 
first  equation  of  (11-  39)  of  TNP  I  gives  for  a  frequency  o  : 


p  =  i  e 

hr  a 


(10  -  11) 


By  eqs.  (9-23)  and  (9-26)  of  the  present  report  we  have 


_  _  _  _  u  iot 

pR  -  m  =  -  *  ee 


(10  -  12) 


Then  it  follows  from  (10-9)  and  (10-11)  that 


n  -  fl  _  _  n  .iot 

Z  0+ft  R  0+fi  £ 


(10  -  13) 


Then  (10-8)  gives  the  nutation  of  the  rotation  axis 


n  =  p  +  n  =  — —  p„ 
°  Kd  -  KR 


( 10  -  14) 
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n  *  -jrree' ot 
R  o+Q 


(10  -  15) 


The  obvious  relation 


p  =  0 
z 


(10  -  16) 


following  from  (10-5)  concludes  the  formulas  for  precession 
and  nutation  of  the  rotation  axis  and  the  z-axis. 

For  the  instantaneous  figure  axis  F  we  have  by  TNP  I, 
eq.  (3  -  27): 


PF  -  f 


c 


(10-  17) 


We  express  c  *  c^  +  ic23  by  (9-2),  introducing  the  secular 
Love  number  defined  by  (9-8).  This  gives 


_  k  u-w 

PF  =  IT  —  ■ 

s 

By  (9-26)  and  (9-28)  this  becomes 


(10  -  18) 


V  x  f  <‘*0  e 


(10  -  19) 


The  nutation  of  the  axis  F  then  follows  from  (10-8)  and 
(10  -  13): 


+  n  = 
z 


(10  -  20) 


We  finally  state  the  relation  which  connects  the  nuta¬ 
tion  number  n  (for  any  axis)  with  the  usual  expression  in 
terms  of  the  Euler  angles  ^  and  6  ,  the  latter  being  the 

obliquity  of  the  ecliptic: 


ae  +  lysine  =  -ine 


(10-21) 


This  is  eq.(ll-55)  of  TNP  I.  The  factor  -i  expresses  a 
rotation  by  -90°  and  has  no  deeper  significance  since  it 
characterizes  only  the  choice  of  coordinate  axes;  the  factor 
e!nt  expresses  the  uniform  rotation  of  the  x°x°x°  system 
to  which  n  refers,  with  respect  to  the  inertial  system  to 
which  $  and  9  refer. 

Forced  polar  motion.  For  the  lunisolar  torque  we  take 
the  expression  (11-40)  of  TNP  I: 


L  -  (C-A)a^  B  e':(wjt  +  Bj) 
j  J 


(10-22) 


as  the  sum  of  contributions 


L.  =  ( C-A)^2B  ,e" 1 j  j ) 


(10-23) 
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of  different  frequencies  uj  with  phase  a .  .  We  put,  as  usual, 
a  =  -Wj  ,  (10-24) 

then 

g-iUjt+gj)  =  giot-igj  ,  e-*6jeiat  >  (10-25) 


So  far  we  have  disregarded  the  phase  factor  e"'6j  because, 
being  common  also  to  u,  v,  w,  it  cancels  in  relations  such 
as  (9-20)  or  (9-42).  We  may  also  take  care  of  this  phase 
factor  by  allowing  uq,  vq,  and  wq  in  (9-3)  and  (9-4)  to 
by  complex  constants.  We  shall  now  use  the  subscripts  j 
to  indicate  that  we  deal  with  contributions  of  the  frequency 

u)  *  • 

J 

We  thus  can  write  (9-28)  as 

w.  =  s2T.e“i(u,jt+6j^  ,  (10  -  26) 

J  J 

whence  ( 9  -  9 )  gi ves 

L.  =  -i ( C-A)fiW . 

J  J 

=  -i(C-A)n2Tje'i(wjt+6j)  • 


(10-27) 


The  comparison  with  (10-23)  shows  that 

Tj  =  i B  j  .  (10  -  28) 

This  relation  is  not  more  mysterious  than  a  rotation  by  90° 
since  iz  is  a  complex  number  obtained  by  rotating  z  by 
the  angle  of  90°  : 

if 

u  =  e  lz  .  (10-29) 

It  expresses  the  fact  that  w  and  hence,  by  (9-42),  u  and 


v  are  normal  to  L  .  Since,  as  we  have  seen  above,  the  polar 
motions  and  nutations,  for  a  given  frequency,  are  all  propor¬ 
tional  to  u  and  w  ,  all  quantities  pR,  p  ,  n_  ,  n  ,and  n 

*  r  R  z  F 


v 

FIGURE  10.2  Forced  nutation  and  polar  motion  for  a  fixed 


frequency 


have  the  same  directions  as  u,  v,  w,  namely  normal  to  L  , 
as  Fig.  10.2  schematically  shows.  Thus  the  points  x° ,  x^ , 

R  and  F  (  the  end  points  of  the  respective  spatial 
unit  vectors  e  ,  e  ,  e  ,  and  e_)  lie  on  a  straight  line. 
Following  (9-26)  and  (9-27),  we  similarly  write 


Uj  -  0£je 


-  i  (to  j  t+8  j  ) 


V.  -  nc.«*i(“Jt+8j)  . 

J  J 
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Then  ( 9  -  29 )  gives 
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(10  -  32) 


note  that  we  have  put  k  =  k.  since  the  elastic  Love  number 
k  also  depends  on  frequency. 

Now  (10-12)  (10-16),  and  (10-19)  give  immediately,  after 

summing  the  contributions  of  the  individual  frequencies: 
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P  forced.  e~i  (ajjt  +  Sj  ) 


p  forced  =  0 
z 


(10  -  33) 


forced  3  l  1J.(e.-T.)e“i(“jt  +  BJ  ) 
j  s  J  J 


Here,  t.  has  the  same  direction  as  t.  =  iB.  . 
J  J  J 

Forced  nutation.  With 


Aw.  =  a.-  =  -(o+n) 

j  j 

eqs  .  (10-13),  (10  -  15)  ,  and  (10  -  20)  ,  combined  with  (10-21), 
gi  ve 


(a  9+iAi/Sin9  )Rforced  =  -i  l  e  .e"  1  ( Au  j t  +  e  j  >  , 

j  J  J 


(a  e+i Arsine  ) 


forced 


■r  Q  -i(Aw.tf8.) 

IwT  Eje  '  1  J  • 

J  J 


(10-34) 


(19tl4*$1n9  )Ffor“d  *  -ll  r^Uj-Tj)  ♦  j2_  c. 


,  -  i  (  Aoj  .  t  +  0  .  ) 


J  J 


For  a  purely  elastic  earth  (no  liquid  core)  we  have 
(9-34)  which  in  view  of  (9-37)  may  also  be  written 
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elastic 

i 


C-A 

nr 


K  .  T  . 

J  J 


C-A 

~r 


1C  .  B  . 
J  J 


(10-35) 


using  (10-28)  and  the  abbreviation  <.  according  to  TNP  I, 
eq.  (6-8).  The  substitution  of  (10-35)  and  of  =  iB^ 
then  immediately  gives  the  forced  (lunisolar)  parts  of  eqs. 
(6-9),  (6-10),  (8-7),  and  (8-11)  of  TNP  I,  which  have  been 
obtained  there  in  a  considerably  more  complicated  way. 

Free  motion.  There  are  two  proper  frequencies,  o1 
correspondi ng  to  the  Chandler  period,  and  o2  responsible  for 
the  "nearly  diurnal  free  wobble"  (NDFW);  cf.  TNP  II,  pp.123- 
126.  Thus  PRfree  have  the  form  of  a  linear  combination 

of  both  frequencies: 

p/ree  =  re10'1  +  m.e1^1  ;  (10-36) 

K  I  L 


a  phase  factor  of  the  form  e1  can  be  incorporated  in  the 
complex  constants  m1  and  .  Then  (10-16)  and  (10-19) 
gi  ve 


n  free  =  0 


k ,  .  k 

r-  m.e  l  +  r~  m  e  2 


free 


(10-37) 


and  being  the  Love  numbers  k  for  the  frequencies 

0,  and  a2  ;  there  is  t  *  0  for  free  motion. 

Free  nutation  is  given  by  (10-13),  (10-15),  and  (10-20) 
with  (10  -  21): 


(Ae+iAtfiSine ) 


free 


«el<Vn,t  ♦ 


Oj+n  1 


+  -1-  m ,e* 
o2+«  2 


(AO  +  iAipsine  )  free  =  i  j  — 1 m  e 1  ^ CT i  +n ^ 1  + 
z  1 


♦  -iL-  m,ei(°2+«>t 

o2+fi  2 


(10  -  38) 


(A0+i Arsine) 


free 


=  -l 


fk, 

IT 


k2 

ir 


a  1  +  S2 


m.e5(o,+fl)t  + 


v“J  ”’e  2  j 


The  nearly  diurnal  free  wobble  (NDFW)  does  not,  so  far, 
seem  to  have  been  confirmed  by  observation,  neither  in  the 
polar  motion  (amplitude  !m  [  )  nor  in  nutation.  For  nutations, 
formulas  (10-38)  show  that  this  amplitude  would  even  be  in¬ 
creased  by  a  factor  on  the  order  of 


1 

►  - 
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a 

a2+a 


-I.0V2V1!  s  500  * 


according  to  TNP  II,  p . 16 7 1 

Thus  we  can  safely  disregard  m2  ,  and  consider  only 
the  Chandler  frequency  .  Then  these  formulas  for  free 

polar  motion  and  nutation  reduce  to  the  free  terms  in  the 
correspondi ng  expressions  for  an  elastic  earth, TNP  I,  eqs.  (6-9), 
(6-10),  (8-7),  and  (8-11),  with  -  aQ  and  k  being  the 
value  kj  for  the  Chandler  frequency. 

Also  the  expressions  (10-33)  and  (10-34)  are  formally 
the  same  as  for  a  purely  elastic  earth,  liquid-core  effects 
entering  only  indirectly:  weakly  through  the  frequency-depend¬ 
ent  Love  numbers  and  more  strongly  through  the  which 

depend  on  core  rotation  5^.  through  (10-32).  This  is  not 
surprising  s.nce  the  basic  relation  for  these  "geometric"  or 
"kinematic"  axes  depend  only  on  the  surface  form  of  our  earth. 

Thus  the  qualitative  picture  of  polar  motion  for  an  elas¬ 
tic  earth  with  a  liquid  core  remains  basically  the  same  as 
for  a  purely  elastic  earth;  cf.  Fig.  6.1  on  p.38  of  TNP  I.  Only 
the  amplitudes  of  forced  polar  motion  are  somewhat  affected 
by  the  liquid  core.  The  orders  of  magnitude  remain  the  same: 

6  m  for  free  polar  motion,  60  cm  for  forced  motion  of  the  ro¬ 
tation  axis,  and  60  m  (!)  for  the  forced  motion  of  the  instanta¬ 
neous  figure  axis;  cf.  TNP  I,  pp.  39-40. 

Finally  we  note  that,  as  far  as  free  polar  motion  is 
concerned,  the  elastic  earth  model,  even  with  a  liquid  core, 
gives  a  highly  schematic  and  unrealistic  picture.  The  actual 
free  polar  motion  is  so  irregular  that  it  can  only  be  deter- 
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mined  by  actual  observations.  On  the  other  hand,  the  forced 
polar  motion  and  nutation  is  represented  well  by  the  elastic 
mantle-liquid  core  model.  Thus  the  formulas  (10-33)and  (10-34) 
describe  them  with  sufficient  accuracy. 
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11.  The  Dynamic  Axes 

These  axes  comprise  the  angular  momentum  axis  and  the 
related  "Celestial  Ephemeris  Pole". 

Angular  momentum  axis.  It  represents  the  direction  of 
the  angular  momentum  vector  H 

Since  the  relation  between  angular  momentum  and  torque, 

dH 

k  "  HI  .  (11-D 

is  independent  of  the  internal  structure  of  the  body  (cf.  TNP  I, 
pp.  48  and  128),  we  have  for  the  nutation  nH  the  same  form 

nH  Cfi(o+aj  (11  '  2) 


as  for  a  rigid  body  or  for  the  Poincarfe  model;  this  is  eq . 

( 13  -  37)  of  TNP  I . 

The  free  nutation  of  the  H  axis  is  zero  since,  in  absence 
of  external  forces,  the  vector  H  remains  unchanged  in  iner¬ 
tial  space.  For  the  forced  motic :  we  have,  summing  over  the 
different  frequences. 


(11  -  3) 


or  by  (9-11),  (10-21),  and  (10-27): 
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Ue+u<|»sine)H*-i£  j  JLL  T  e"'  ( t+6j  >  .  (11-4) 

j  j  J 

Using  (10-28)  we  have  the  alternative  form 


(Au  +  iAijiSine  )  H  =  £  £  B  e-,{Awjt+6j^  . 

j  j 


(11-5) 


This  is  identical  to  the  nutational  part  of  eq.(7-14)  of 
TNP  I,  but  the  derivation  is  much  simpler  in  the  present  way. 

It  is  also  identical  to  the  nutation  of  H  as  given  by  TNP  I, 
eq.(13-69)  for  the  Poincarfe  model. 

Polar  motion  for  the  H  axis  may  be  obtained  from  (10-8): 


PH 


(11  -  6) 


Using  (10-13)  and  (11-4)  and  summing  over  different  frequen¬ 
cies  we  get 


PH 


fA  JLl 

A  w  j 


T  . 

J 


cl  e'i(ujt+V 

Ji 


(11-7) 


A  thorough  check  is  provided  by  computinq  p..  from  an 
expression  of  the  vector  H  in  the  body  frame.  The  angular 
momentum  equation 


k  V''  V* 


108 


H  =  C  u,+  h  (11-8) 

(TNP  I,  p . 1 3 )  holds  for  an  arbitrary  nonrigid  body  (even  with 
a  liquid  core),  provided  the  relative  angular  momentum  ji  is 
taken  into  account  besides  the  main  term  U  ,  C  being  the 
inertia  tensor  and  w  the  rotation  vector  as  usual. 

For  Tisserand  axes  referred  to  the  mantle,  the  relative 
angular  momentum  h  is  due  to  core  rotation.  Thus 


(H-9) 


since  Xj=  0  by  (8-9);  we  have  disregarded  the  second-order 
products  of  the  anomalous  tensor  of  inertia  for  the  core. 


cc 


with  v.  as  usual 

A  i 


The  linearization  of  Zu  is  provided  by  eq.(3-22)  of 
TNP  I;  by  (11-8)  we  have  to  add  h_  as  given  by  (11-9).  The 
result  is 


H  =  Afim  +  nc  +  Ac(Xl  +  ix2)  • 
By  ( 1  -  3 )  this  is 

H  =  Aiim  +  ac  +  A  v  . 


(11  -  10) 
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Here 


H  *  H1  +  iH2  (11-11) 

is  the  complex  number  combining  the  x1  and  x2  components 
of  the  angular  momentum  vector  H  .  The  components  of  the 
correspond!- ng  unit  vector  are  obtained  by  dividing  by  its 
(approximate)  length  Cn  ,  whence 


PH 


(11  -  12) 


we  obtain  indeed  p  since  everything  is  referred  to  the  body 

H 

frame  X1X2X3  •  witl1  (9  -  23)  and  (10-17)  this  may  be  written 


After  substituting  p  from  (10-12),  p  from  (10-19)  and 

R  r 

v  from  (9-20),  and  neglecting  second-order  terms  we  get 
(11-7),  which  completes  the  check. 

This  equation  can  also  be  used  to  compute  the  free  polar 
moti on  of  H  .  Core  rotation  only  plays  a  role  in  the  tidal 
frequency  range,  that  is,  for  .  =  a  ;  its  effect  is  rather 
small  for  the  Chandler  frequency,  cf.TNP  I,  p . 1 3 3 .  Thus,  since 
the  calculation  of  free  motion  is  of  little  practical  value 
anyway,  we  may  neglect  v  in  (11-13),  obtaining 


„  free 
PH 


A  free 
Z  PR 


+  ZJl  P(:free 


(11  -  14) 


Using  (10-36)  and  (10-37)  and  neglecting  the  hardly  observ 
able  NDFW  (m  =  0)  ,  we  get 


A  .  C-A  ki 

IT  +  ~T~  r~ 


_  -ic.t 
m  j  e  l 


(11-15) 


the  same  as  for  a  purely  elastic  earth  (TNP  I,  p.37). 

The  Celestial  Ephemeris  Pole.  The  corresponding  axis 
has  been  adopted  at  the  General  Assembly  of  the  International 
Astronomical  Union  in  Montreal  in  1979  to  define  the  celestial 
pole  for  reference  purposes,  according  to  Commission  Resolu¬ 
tion  ( 3 ) : 

"Commissions  4,  8,  19  and  31 

endorse  the  recommendations  given  in  the  Report  of  the 
Working  Group  on  Nutation,  as  set  out  below,  and  recom¬ 
mend  that  they  shall  be  used  in  the  national  and  inter¬ 
national  ephemerides  for  the  years  1984  onwards,  and  in 
all  other  relevant  astronomical  work. 

Recommendations  of  the  Working  Group  on  Nutation 

Whereas  the  complete  theory  of  the  general  nutational 
motion  of  the  Earth  about  its  centre  of  mass  may  be  de¬ 
scribed  as  the  sum  of  two  components,  (i)  astronomical 
nutation,  commonly  referred  to  as  nutation,  which  is 
motion  with  respect  to  a  space-fixed  coordinate  system, 
and  (ii)  polar  motion,  which  is  motion  with  respect  to 
a  body-fixed  coordinate  system,  it  is  recommended  that: 
(a)  astronomical  nutation  be  computed  for  the  "Celestial 
Epnemeris  Pole"  using  a  non-rigid  model  of  the  Earth 
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such  that  there  are  no  nearly  diurnal  motions  of  this  celes¬ 
tial  pole  with  respect  to  ei ther  space-fixed  or  body-fixed 
coordinates,  which  can  be  calculated  from  torques  external  to 
the  Earth  and  its  atmosphere  , 

(b)  the  numerical  values  given  in  Table  1  of  the  complete  re¬ 
port  be  used  for  computing  astronomical  nutation  of  the  "Celes¬ 
tial  EDhe men's  Pole".  " 

Cf.  (IAU,  1980,  pp.  40-41).  For  the  scientific  background  cf. 
(Leick  and  Mueller,  1979;  Moritz,  1979;  Mueller,  1981)  and  also 
TNP  I,  pp.  53-59. 

The  Celestial  Ephemeris  Pole  C  thus  corresponds  to  the 
angular  momentum  axis  freed  from  the  lunisolar  diurnal  motion, 
seen  from  the  earth-fixed  body  frame,  that  is  with  respect  to 
polar  motion.  With  respect  to  the  inertial  frame,  it  shares 
with  H  the  property  that  its  free  nutation  (which  has  a 
nearly-di urnal  period)  is  zero.  Thus  in  fact,  C  contains  no 
diurnal  motions,  neither  in  the  body  frame  nor  in  the  inertial 
frame,  as  the  above  resolution  requires. 

Hence  the  free  polar  motion  of  the  Celestial  Ephemeris 
Pole  C  equals  the  free  part  of  p  : 

Pc  =  PHfree  .  (11  -  16) 

An  analytical  expression  is  (11-  15),  which  is  of  little  prac¬ 
tical  use,  however,  since  free  polar  motion  is  too  irregular 
to  be  analytically  predicted.  So  pc  must  come  from  obser¬ 
vation. 

As  far  as  the  nutation  is  concerned  we  have  seen  above 
that  its  nearly-di urnal  part  is  zero.  There  remains  the  luni¬ 
solar  analytical  part  which  can  be  analytically  described  and 
predicted  well.  Free  motion  playing  no  part,  the  pole  C 
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does  not  depend  on  polar  motion.  Its  nutation  is  therefore  the 
same  whether  polar  motion  exists  or  is  zero.  In  the  latter  case 
C  coincides  with  the  z-axis.  Cf.  TNP  I,  Fig.  6.1  on  p.38: 
if  there  is  no  polar  motion,  then  C  »  H  coincides  with  the 

O 

origin  0  representing  the  z-axis.  See  also  TNP  I,  pp.  59-60. 

Thus  the  nutation  of  C  equals  the  forced  (lunisolar) 
nutation  of  the  z-  axis  as  given  by  (10-34): 

Ufl  +  U*)c--iJ  e"!(A“.it+Bi)  (11-17) 

j  j 

Concluding  Remark.  Only  the  forced  (lunisolar)  parts  of 
polar  motion  and  nutation  can  be  accurately  described  by  the 
elastic  mantle-liquid  core  model  and  hence  accurately  predicted 
Therefore,  the  main  practical  result  of  the  present  report  are 
the  formulas  (  10  -  33),  (  10  -  34),  (  1 1  -  5  ) ,  ( 11-7) ,  (11-17). They 
contain  the  coefficients  e.  and  x.  which  are  given  by 
(  10  -  23),  (  10  -  28),  and  (10-32). 
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